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A USERS MANUAL 

FOR COMPUTERIZATION 

OF HYDROGEOLOGIC DATA 

FOR A RIVER ALLUVIUM 

By James \V. Nflney, Bill B. Barnes, and Douglas C. Kent 1 

ABSTRACT 

A system for rapid storage and retrieval of data from a multilayered 
aquifer based on the hydraulic properties of the aquifer was developed. 
The manual consists of the data-retrieval system, control- card setups and 
specifications, lists of variables, formats and source listings for data stor- 
age and retrieval, and graphical illustrations demonstrating how the system 
works. Data for the examples were collected in the alluvium of the 
\Vaslnla River near Chickasha, Okla, Those data were processed according- 
to the manual to produce a lithologic cross section, an isopachous map, 
a specific snbdatum map, a distributed test-hole-data map, and a perme- 
ability-distribution map. The system proved capable of storing 1 and re- 
trieving- a large volume of layered hydrogeolog'ic data, and it will have 
wide geographic application where hydrogeolog'ic studies are undertaken. 
The system may also be adapted to other types of materials because it 
uses numerical coding' exclusively to describe the hydrog'eolog'ic data to be 
stored, retrieved, and manipulated. KEYWORDS: computerization of 
hydroffcologic data, ground water, ground -water storage, ground-water 
transmissibility, hydrogeology, hydrology, lithology, soil permeability. 



INTRODUCTION 

A layered lithologic system created the need 
for specialized computer programs to store, 
retrieve, and manipulate hydvogeologic data. 
Existing programs (7, 5. 7) J would store data 
on labeled punched cards or provide for the 
alphameric coding' of data on punched cards 
and store and retrieve ground-water hydro- 
graphs. This manual, however, has been 
developed for the selection and retrieval of 
ground-water data based on the hydraulic prop- 
erties of an alluvial aquifer. It documents tech- 
niques described in an earlier paper by Kent et 



1 Geologist and computer specialist, Agricultural Re- 
search Service, U.S 1 . Department of Agriculture, 
Chickasha, Okla. 73018, and associate professor of 
geology, Oklahoma State University, Stilhvater, Okln. 
74074. 

- Italic numbers in parentheses refer to items in "Lit- 
erature Cited" preceding appendix A. 



al. (.j*) ; excerpts from that paper are included 
here for convenience. 

The retrieval system was applied to define 
boundaries and internal hydraulic character- 
istics of an aquifer in the Waslrila River al- 
luvium near Chickasha, Okla. (fig. 1). The 
alluvium is characterized by several discon- 
tinuous layers of silty clay, sand, and gravel 
similar to those found in stratified alluvial fans, 
as well as some stratified basin and Coastal 
Plain sediments. Although the approach used 
in this study was specifically designed for al- 
luvial aquifer systems, as shown by Naney and 
Kent (fi), it can be used for other similar 
aquifers. 

The computer codes presented in this manual 
are written in FORTRAN IV for the IBM 360- 
65 computer. Some of the techniques used for 
processing the data in the examples were de- 
veloped on the IBM 1130. A restriction on the 
word length for coding some of the data was 
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FtfiURE l. Hydiiiulic-coefficient ranges in the Washita River alluvium. 



encountered with thy IBM 1100. The adaption 
of the>e programs to both computers, however, 
lias expanded the availability of the technique 
for which tins manual is written. 

Three types of data are documented and 
stored in this data-handling system; (1) gen- 
eral test-hole information {test-hole data); 
(2) stratification characteristics (layer data); 
and (o) water-level record^ (water-level 
data). The test-hole data consist of the well 
number, location, elevation of ground surface, 
elevation of top of pipe (uell casing) , the num- 
ber of layers and the method of drilling. Layer 
data include layer number and thickness, meth- 
od of analy.sis, lithologic type and color, and 



approximate range of the hydraulic coefficient 
for each layer. The water-level data include well 
number, record type, date, time, and water-level 
elevation. 

DATA-RETRIEVAL SYSTEM 

The data-retrieval system (fig. 2) is based 
on descriptions of the lithology in terms of the 
type and color of material and of the hydraulic 
properties associated with each layer. These 
descriptions are stored by test-hole number, 
latitude and longitude, and layer number on 
a disk or magnetic tape. General test-hole in- 
formation and water-level records are stored 
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FIGURE 3. Flow chart for retrieval of hydrogeologic data. 



by well number and location on separate disks 
or magnetic tapes. Specific information related 
to each of these data is numerically coded and 
stored by variable number. Combinations of 
these stored data are used as criteria in the 
selection of other data or as data that are to 
be shown graphically. 

The hydraulic characteristics of a multi- 
layered aquifer system can be evaluated by 
options that allow selective retrieval of data. 
Ten options for data selection and comparison 
are specified on control cards. A general flow 
chart of the retrieval program is shown in 
figure 3. Figure 4 shows how each specific 
option is selected. Control-card layouts for 
retrieval options A through J and a description 
of card types A through N appear in appendix 
A. The variable name and columns assigned 
to each code are shown for each option. A list of 
variable names with descriptions also appears 
in appendix A. In addition to the options de- 
scribed, subroutine packages may be introduced 
to be called by specified control cards. 

Data-storage format. The data-storage pro- 
gram reads data from punched cards in the 
formats shown in appendix B and stores them 
on magnetic tape. Twenty-one variables are 
stored for each test hole, including variables 
that describe the entire well, such as the type 
of well casing or the elevation of the top of 
the pipe. These variables are presented in table 
1. 

Twenty lithologic zones or layers are allowed 
for each test hole. Currently, no more than 14 
distinguishable layers have been stored for 
any one test hole. There are eight variables 
describing the litholog'ic characteristics of each 
layer, such as layer thickness or the type of 
material. These variables (table 1) are read 
from punched cards according to the format 
shown in appendix B and stored on magnetic 
tape. 

Ground-water levels are stored on magnetic 
tape for use in conjunction with the test-hole 
and lithologic-layer data. The format used on 
punched cards to store the data is shown in 
appendix B, Storage and retrieval programs 
for test-hole, layer, and water-level data, along 
with comments and execution cards, appear in 
appendix C. The three basic graphical forms 
produced are test-hole cross sections, lithology- 
distribution maps, and isopachous maps. A 
separate retrieval program (not shown in fig- 



ure 2) provides data for maps that show the 
distribution and change in level of the ground- 
water surface. 

Coding and storage of lithohffic data. The 
lithologic type and color are described with 
unique combinations of numbers (table 2). Only 
five numbers representing grain size are used in 
the code to describe the lithologic type because 

TABLE 1. -Data-processing code for tape stor- 
age, entire well 



Variable 



Term 



Code form 1 



1 Well number or test-hole 

number XXXX 

2 Watershed number XXX 

3 Location (latitude; degrees, 

minutes, seconds) XX XX XX 

4 Location (longitude; degrees, 

minutes, seconds) XX XX XX 

5 Top of pipe elevation XXXX.XX 

6 Ground-surface elevation XXXX.XX 

7 Bedrock elevation XXXX.XX 

8 Total depth of test-hole log (ft) XXX 

9 Distance from logged well to 

observation well (ft) XXXXXX 

10 Total depth of observation well 

(ft) XXX 

11 Drilling method X 

Core 1 

Rotary 2 

Split spoon 3 

Auger 4 

12 Casing type X 

Metal 1 

Plastic 2 

18 Casing length (ft) XXX 

14 Casing diameter (in) XX.XX 

15 Screen XXX.XX 

Screen length (ft) XX.XX 

Screen diameter (in) . . XXX.XX 

16 Pump type X 

Submersible 1 

Jet 2 

Turbine 3 

Centrifugal 4 

17 Pump use X 

Domestic 1 

Irrigation 2 

Experimental 3 

18 Pump capacity (gal/min) - XXXX 

19 Number of stages XX 

20 Pump diameter (in) XX.X 

21 Number of layers XX 

1 An X represents a space provided for the variable- 
number code. A number indicates a specific descriptive 
adjective to bo coded in the space provided for that 
variable, For example, a 4 put in the space designated 
by the X after "drilling method" would only select wells 
drilled by an auger. 
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FIGURE 4. Flow chart for 

logs, samples are taken from holes cored ad- 
jacent to formerly drilled test sites. Grain size 
of each core sample is visually described and 
analyzed for comparison with descriptions 
based on drill cuttings. The median grain size 
is equivalent to the grain size having the high- 
est percentage of occurrence in the drill-cutting 
sample and is termed "primary constituent." 
A gram size representing a smaller percentage 
of the sample is labeled "secondary constitu- 
ent." With these comparisons and terms, quanti- 
tative relationships can be transposed into 
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retrieval control cards, A-N. 

terms for coding the visual description of drill 
cuttings. 

Hydraulic properties of the aquifer are de- 
termined by using a relationship between grain 
size and the hydraulic coefficients (fig 1 . 6). 
The hydraulic coefficients (the coefficients of 
permeability and storage) are based on field- 
pump tests of selected alluvial material, lab- 
oratory tests of the core samples, and grain 
size-specific yield relationships (2). The hy- 
draulic properties are coded by dividing the 
primary and secondary grain sizes into four 



identifiable ranges (table 3). Undifferentiated 
sand is included in a fifth range. 

A computer program (appendix D) is used 
to assign a range number (3-5) to the layers 
of each well on the basis of the coded lithologic 
description shown in table 4. Hydraulic-coef- 
ficient values representing these ranges can be 
arbitrarily selected from a graphical relation- 
ship similar to the one shown in figure 5. The 
grouping technique makes them applicable to 
all drill-log 1 data, with only the values repre- 
sented by each code being subject to change 



TABLE 2, Data-processing code for tape storage, stmtigraphic layers 



Variable 



Term 



Code form 1 



Variable 



Term 



Code form 1 



1 Layer number XX 

2 Layer thickness -- ... XX 

3 Type of material:- 1 column 1, num- XXXXX 

her of primary constituents (zero 
if unknown) ; column ~, code of 
prmiaiy constituent; column 3, 
code of adjective for primary 
constituent (zero if primary 
constituent has no adjective) ; 
column !,, numeral 1 if second 
primary constituent or secondary 
constituent has no adjective 
(zero if no secondary constitu- 
ent) , or code of secondary con- 
stituent if it has an adjective, or 
code of second primary constitu- 
ent with adjective; column 5, 
code of second primary constitu- 
ent with no adjective or second- 
ary constituent with no adjective 
(zero if no second primary con- 
stituent oi 1 secondary constituent 
exists), or code of adjective for 
secondary constituent, or code of 
adjective foi second primary 
constituent. 

Constituents X 

Clay - . , . ... 

Silt 1 

Sam! 2 

Gravel 3 

Sandstone 4 

Carbonate ... . , . . 5 

Shale 6 

Evapoi-ite - - 7 

Untlifferentiated alluvium .- 8 

Undifferentiated bedrock 9 



3 Type of material Continued ;- 

Log-description adjective X 

Very fine 2 

Fine 3 

Medium 4 

Coarse 6 

Lime fi 

Dolomite 7 

Halite 8 

Fossiliferous 9 

4 Color of material: major color is XX 

the left digit, minor color the 
right digit (a zero in either place 
indicates no color).'' 

5 Analytical method X 

Visual accumulation tube 1 

Sieve 2 

Driller inspection (field) 3 

Microscopic examination . . 4 
HYDHAULIC COEFFICIENTS! 

6 Transmissibility (T) (gpd/ft) X 

Values from to 8,000 ... - 1 

Values from 8,000 to 30,000 . 2 

Values from 30,000 to 80,000 . 3 

Values from 80,000 to 160,000 4 

7 Permeability (K) <gpd/ft=) X 

Values from to 80 . . 1 

Values from 80 to 300 - - . . 2 

Values from 300 to 800 ... 3 

Values from 800 to 1,500 ... 4 

8 Storage coefficient (S) (dimension- 

less ratio) X 

Values from 0,0000 to 0.0001 1 

Values from 0.0001 to 0.0009 . . 2 

Values from 0.0009 to 0.0020 . 3 

Values from 0.0020 to 0.0040 - 4 



1 An X represents a space provided for the variable-number code. A number indicates a specific descriptive ad- 
jective to be coded in the space provided for that variable, 

- For example, a log desciiption of a primary and a secondary constituent as "medium sand, with some fine sand" 
results in a code name of "sand, medium, with some sand, fine," or a number code of 12423. Again, a log- descrip- 
tion of two primary constituents as "medium sand, and fine sand" results in a code name of "sand, medium, and 
sand, fine," or a number code of 22423. 

1 For example, a sample with a major color of red biit no minor color would be coded 20, one with no major 
color but a minor color of red, 02; a sample with a major color of red and a minor color of brown would be coded 
25, one with a major color of brown and a minor color of red, 52. A code of 00 indicates no color at all was re- 
corded. 

1 Code numbers 1-4 represent a range of values for T, K, and S, 



based upon the geological materials being- 
investigated. 

Retrieval o/ litlwlofjic data. Data can be re- 
trieved in many forms, but the most useful are 
cross sections and maps that define aquifer 
boundaries and the horizontal or vertical dis- 
tribution of the hydraulic properties in the 
aquifer. The schematic diagram of the retrieval 
system (fig. 6) shows the assignment and use 
of the ranges and corresponding estimated 
values of the hydraulic coefficients for presen- 
tation m different graphical forms. Retrieval 
and testing of test-hole, layer, and water-level 
data are essential to provide information neces- 
sary for selection and plotting of specified 
ranges of the hydraulic coefficients (fig. 2). 
Retrieval of test-hole data provides information 
necessary for selection and plotting of the 
ranges at a specified depth or subdatum eleva- 
tion. Water-level data are used to select ranges 
characteristic of each layer within the satu- 
rated zone. Testing on the range number of the 
hydraulic coefficient permits the selection of 
specified ranges of the hydraulic coefficients. 



TABLE 3. Litholoyic type and associated mnye 
of the hydraulic coefficients 



Range 1 



Lithologic description of log j 



1 Silt;- silt with some very fine sand; veiy 

fine sand with some silt; silt ami very 
fine sand. 

2 Very fine sand; very fine sand with some 

fine sand; fine sand with some very 
fine sand; very fine .sand and fine sand. 

3 Fine sand; fine sand with some medium 

sand; medium sand with some fine 
&and; fine and medium sand. 

4 Medium sand; coarse sand and gravel; 

medium sand with some coarse sand and 
gravel; coarse &mid and gravel with 
some medium sand; medium and coarse 
sand. 

5 Undifferentiated material in ranges 2, 3, 

and 4, 

1 Zero is used when a range cannot be as&iffned to a 
layer. 

- Combination of primary and secondary constituents 
that are classified as indicated. 

1 Silt and clay arc undifferentiated, 



TABLE 4. List of codes representing four groups of relative permeability 



Sand 


Clay and silt Very fine to . 
/I) f Fine to medium 

(2) (3 > 


Medium to 
coarse 
(4) 


10000 20011 12211 22211 12324 22324 


12500 22524 


10011 20012 12311 22311 12423 22423 


12524 22425 


10012 20013 12011 21122 12300 


12426 23125 


1001,1 11122 21123 12400 


13000 23124 


20017 11123 22223 12413 


13012 22513 


10017 2001!) 12200 22322 


13125 


1001!) 23010 12300 


13124 


10125 21010 12223 


12513 


10124 21012 12322 




10123 21013 




10122 




12210 




12310 




12410 




12510 




12010 




13010 




11000 




11010 




11012 




11013 




11017 




11019 




13011 




1GOOO 





10,000 




FiGUHE 5. Relationship between hydiauhc-coefficient ranges, median grain si?,e, and 

the coefficient of permeability. 



One important aspect of the retrieval system 
is the use of water-level records for selecting 
data in the saturated zone. A subroutine which 
computes a weighted average of all measure- 
ments within a selected time period is used to 
determine the average elevation of the ground- 
water surface in each well. 

GRAPHICAL PRESENTATIONS 

Examples are given of five basic types of 
graphical presentations of data retrieved. One 



cross-section display and four possible map- 
view displays are shown. A flow chart of the 
steps followed in preparing graphical displays 
and the list of variable names, with descrip- 
tions, used to plot selected data are shown in 
appendix E. 

Cross sections of drill-hole logs. -Drill-hole 
data are punched on cards and may be used 
as input data for plotter programs. The con- 
trol-card setup required to retrieve data for the 
plotting of cross sections is shown in figure 7. 
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FIGURE 0. Data-retrieval system using hydraulic coefficients, 
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FIGURE 7. Control cards for retrieval of data for plotting cross sections. 
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For plotting' options available, a description of 
each retrieval control card used is found in 
appendix A. 

A plotted cross section of drill holes near 
Alex, Okla., is shown in figure 8. Appendix F 
shows the data cards used to plot cross sections. 
The variable name, type of plot, and column 
in which each variable is coded are shown for 
each card. These data may then be used on the 
plotting equipment available to the user. 



Maps of distributed drill-hole data. The 
total depths of the test holes were plotted as a 
distributive map. The control cards for the 
plotter program are listed in figure 9, and a 
map showing the distribution of the total depth 
of each test-hole log appears in figure 10. The 
map shows the depth to which lithologic data 
are available, since some test holes were not 
drilled to bedrock. Appendix F contains the 
data cards used to plot selected data maps, with 
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FIGURE 11. Contiol cards for retrieval of hydraulic-coefficient ranges at a specified 

sub datum. 
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FIGURE 12. Map of hydraulic ranges for a specified suhdatum elevation. 



the type of plot, variable name, and column 
containing the variable code indicated. Any of 
the 21 test-hole variables may be retrieved for 
plotting as a distributive map. Appendix E 
gives a list of variable names and their descrip- 
tions for all example plots. 

Maps of data from specified elevations. 
Because hydrogeologic data concern informa- 
tion that must be gathered at some depth below 
the surface of the ground, a convenient method 
of looking at various levels, or subdatum eleva- 



tions, was developed by the authors. Digital 
computer techniques are used to develop a 
planar surface of lithologic data for a specified 
elevation above mean sea level. Figure 11 shows 
the control cards for the subdatum searching. 
A map of the test holes that have lithology with 
hydraulic coefficients ranging from 1 to 5 at 
the specified subdatum is presented in figure 
12. Appendix F contains the data cards needed 
to plot the map shown in figure 12. 
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FIGURE 13. Control cards for retrieval of weighted averages of the hydtaulic co- 
efficients of all layers in each hole. 
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FIGURE 14, Map of iwevage permeability coefficients. 



Maps of average coefficients of permeability 
and storage, When the ground-water flow 
through a specified cross section or the amount 
of ground water stored in a region of ti water- 
shod is needed, it is sometimes desirable to 
average the coefficients of permeability and 
storage ol the layers that the test hole pene- 
trates. Figure 18 shows the control cards for 
averaging these hydraulic parameters. 

Those averaged values of K (permeability) , 
when displayed in map view (as shown in fig- 



ure 14) , provide a distribution of test holes 
that are potentially high or low in ground-water 
yield. A distributive map of S (storage co- 
efficient) may be produced using the system. 
It is possible to determine suitable well field 
locations and possible artificial recharge areas 
using aveal maps (fig. 1) developed from the 
distributed test-hole data. 

The retrieved data cards used to plot the 
map shown in figure 14 are presented in 
appendix "P. 
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15. Control caidh for retiieval of an ibopachous map of hthology having 
peimealiility ranges 1-5. 
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FifiURE 16, Isopachous map of hythaiilic properties, ranges 1-6. 



Isopachous maps of specific hydraulic prop- 
erties. An isopachous map of all the lithologies 
within an alluvial system having the same 
hydraulic properties was also plotted. Such a 
map is used by planners when well fields are 
being installed or augmented. The amount of 
ground water available in an area may be de- 
termined from such a map (3) . 



A control-card setup for retrieving data for 
an isopachous map of all subsurface materials 
having a hydraulic code 1-5 is presented in 
figure 15. The isopachous map plotted is shown 
in figure 16. 

Appendix F shows the retrieved data cards 
used to plot the map shown in figure 16. 
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APPENDIX A.-CONTROL-CARD SETUPS, CARD SP 

AND LIST OF VARIABLES FOR OPTIONS A-N 



Control-Card Setups 
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FIGURE A-l. Retrieval option A; special drill-hole data. 
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FIGURE A-2. Retrieval option B : type of material. 
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FIGURE A-3, Retrieval option C; color of material. 
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FIGURE A-'}. Retrieval option D : data from specified subdatum planes. 
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FIGURE A-5. Retrieval option E: data from selected heights above bedrock. 
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FIGURE A-G. Retrieval option F: data for logging a test hole. 
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FIGURE A-7. Retrieval option G . data from saturated zone. 
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FIGURE A-8. Retrieval option H: data from specified depths below ground surface. 
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FIGURE A-9. Retrieval option I : average hydraulic coefficients for the saturated zone. 
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FIGURE A-10, Retrieval option J: selection of all test-hole and layer data, 
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FIGURE A-ll, Description of card types A, B, and C. 
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FIGURE A-l 2. Description of card types D, E, and F. 
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FIGURE A-13. Description of card types G, H, and I. 
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FIGURE A-14. Description of card types J, K, and L. 
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FIGURE A-16. Descnption of card types M and N, 
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List of Variables 

Description of variable 
Card Type A 

Number of test holes to be searched for data 
retrieval. 

Number of lithologic-layer variables to be 
retrieved. 

Number of selected values of lithologic-layer 
variables to be retrieved, 

Number of additional lithologic-layer data 
to be retrieved. 

Number of descriptive well data to be re- 
trieved. 

Number of lithologic materials to be re- 
trieved. 

Number of lithologic colors to be retrieved. 

If a 1 is punched in column 40, data are 
retrieved from subdatum planes (see card 
type J). If column 40 is blank, this 
option is bypassed. 
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Vai table 
numa 



Option 



Format 



Description oj variable 



INC 



ILOG 



ILVL 



IGSUB 



KTESK 



KTESS 



KSAT 



IPUNII(l) 

IPUNH<2) 

IPUNH(3) 
IPUNIT(4) 
IPUNH(B) 
IPUNII(6> 

Note: 
IWELL(n) 



G 



H 



D,E,G 



B-J 



C,F 



B 



I 



Caid Type A Continued 

15 An integer punched in columns 41 15 is the 

number of inciements to be searched 
above bedrock (see card type K) . If 
columns 41-45 are blank, this option i^ 
bypassed. 

IB A 1 in column 50 indicates a total test-hole- 

log 1 retrieval. If column 50 is blank, only 
the selected layer variable is retrieved. 

15 A 1 in column 55 calls water level for 

defining saturated zone (see caid type 
Nj. If column 55 is blank, this option is 
bypassed. 

15 An integer in columns 50-00 indicates the 

number of feet below ground bin-face to 
be searched by increments. If columns 
56-60 are blank, this option is bypassed. 

15 An integer in columns 01-05 indicates the 

number of permeability values assigned 
(see card type L). If columns Gl-65 are 
blank, this option is bypassed. 

IB An integer in columns GG-70 indicates the 

number of storage-coefficient values as- 
signed (see card type M). If columns 
GG-70 are blank, this option is bypassed. 

15 If column 75 is a 1, K and S values are 

averaged for the total \vcH log. If column 
75 is blank, K and S values are averaged 
only for the saturated zone. 

Card Typo B 
IB 



A 1 in column 5 causes retrieved layer data 
to be punched on cards. 

A 1 in column 10 causes the net thickness 
of selected lithologic layers to be punched 
on cards. 

A 1 in column 15 causes test-hole (well) 
variables to be punched on cards, 

A 1 in column 20 causes lithologic color to 
be punched on cards. 

A 1 in column 25 causes lithologic material 
type to he punched on cards. 

A 1 in column 30 causes the weighted aver- 
age of hydraulic coefficients to be 
punched on cards. 
A blank in appropriate columns will bypass punch output for that option. 
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15 



IB 



IB 



IB 



IVAR(w) 



Card Type C 

A-J 1615 This array is used to select test-hole data 

by numbers. As many as 1,000 drill-hole 
numbers may be retrieved from this array 
selectively. 
Note: indicates array size. 

Card Type D 

A 815 Data are retrieved from this array for as 

many as 8 lithologic-layer variables. 
These variables describe hydraulic and 
geologic characteristics for discrete 
layers. 
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Variable 
iiame 



TVALUOO 



Option 



Format 

Card Type E 
GP12.4 



Description of variable 



Data are retrieved from this array for as 
many as 8 test values (G per card). One 
or more values may be tested for each 
variable indicated on card type D, Addi- 
tional card type E is needed for more 
than G values. 



IADVR(n) 



C-J 



Card Type F 
815 



Data are retrieved from this array for as 
many as 8 layer-variable values in addi- 
tion to the ones tested for using' card 
types D and E. 



ISTOR(a) 



A-J 



Card Type G 
1615 



Data are retrieved from this array for as 
many as 21 drill-hole variables. An ad- 
ditional card, type G, is needed for more 

than 16 test-hole variables. 



MAJTS(w) 



Card Type H 

1015 Data are retrieved from this array for as 
many as 10 discrete types of material 
based upon size-distribution analyses. 



MACOL(tt) 



Card Type I 

1015 Data are retrieved from this array for as 
many as 10 discrete colors of material. 



SUBDA 



SUBTS 



SUBIN 



Card Type J 

F10.0 Maximum elevation to be searched (feet 
above mean sea level) when selecting data 
from horizontal subdatum planes. 

F10.0 Minimum elevation to be searched (feet 
above mean sea level) when selecting data 
from horizontal subdatum planes. 

F10.0 The distance in feet to be stepped down for 
subsequent subdatum planes from 
SUBDA to SUBTS. 



AINC(H) 



AKVAL(w) 



SVAL() 



Card Type K 
1615 



Data are retrieved from this array for as 
many as 16 increments (feet above bed- 
rock) for which data are retrieved. 



Card Type L 
10F8.0 



Data are retrieved from this array for as 
many as 10 discrete permeability values. 



Card Type M 

10F8.0 Data are retrieved from this array for as 
many as 10 discrete storage-coefficient 
values. 
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Variable , 

namc Option Format Description of variable 

Card Type N 

IMO G,I 15 Month of beginning' of time period for which 

ground-water levels are averaged when 
data are retrieved from the saturated 
zone. 

G,I 15 Day of beginning of time period. 

G,I 15 Year of beginning of time period. 

ITIM G,I 15 Time (military) of beginning of time 

period. 

JMO G,I 15 Month of ending; of time period. 

JDY G,I 15 Day of ending of time period. 

JYR G,I IB Year of ending of time period. 

JTIM G,I 16. Time (military) of ending of time period. 
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APPENDIX B.-FORMATS USED FOR DATA STORAGE 



INPUT -CARD FORMAT USED TO STORE DRILL-HOLE DATA 



r? 


IS 


IS 


I 5 


15 


FII.3 


FM.3 


FII.3 


FII.3 






















COLUMN 67-80 






















i 


















BLANK 





o 

z 


o 

z 


O 

z 


d 

z 


Ul 


Ul 


1 


1 






HI 


ui 


UJ 


H 


S 


> 








UJ 


o 


o 


o 


S 


3 

m 


3 


CD 


3 






5 


E 

Q 


DRILL 


1 
g 


g 


1 


o: 

s 


i 


3 





INPUT-CARD FORMAT USED TO STORE LAYER DATA 



14 


F40 


F4.0 


X, F5.0 


F3.0 


T 3.0 


F7.0 


F7.0 


F 10,0 


14 




a: 


















r 




a 


















a 




i 


















u 




Eg 


















X 


COLUMNS 53 -BO 


*g 


















s| 


BLANK 


UJ 
















_ 


aw 




li 
















I- 


la 












_j 


o 


P 










1 






a 


-J 

< 




i 


5 


e 


* 


u 

E 







z 
u 


i 


z 
* 
o 




UJ 

1 

U. 


u, 
o 


s 


j 
m 


Sg 


e 

3 


o 

Ul 


3 
fe 







* 




O 




Q 


5 


3 


o 


(K 




g 


5 


LAYER 


UJ 
0. 


COLOR 


fi 


I 


i 


3 

V) 


HI 






ID 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 







INPUT -CARD FORMAT USED TO STORE WATER-LEVEL DATA 



14 


12 


17 


13 


F7.2 


17 


15 


F7.2 


17 


15 


F7.2 












{ 




















a. 




DC 














COLUMN 64-60 






5g 


W 
















BLANK 




4 


S 


9 


g fc u 
















* 


* 


>-* " 


P 


* Q 
















d 


d 


2 


j. 


S u "" 
















z 


Z 


t 


tc 


fp sk 






























X 










-j 

Ul 


Ul 

o 
o 
u 


s- 


13 

5 


Q U, R 


i 


? 


t 

Ul 


S 

a 


Ul 

F 


i 





ALL I FORMAT INPUT DATA IS RIGHT JUSTIFIED 
* CODE I RECORDER-CHART DATA CODE 2 -TAPE- DOWN DATA 
* NO ZERO PUNCHES NEEDED IN BLANKS ON CARDS 



FIGURE B-l. Input card formats. 
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APPENDIX C. 



-SOURCE LISTING FOR STORAGE AND RETRIEVAL 
PROGRAMS 



CARD 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 



42 
43 
44 
45 
46 
47 
48 
49 
50 

51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

99997 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

99998 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 



APPENDIX C 



RETRIEVAL PROGRAM FOR HYDROGEOLUGIC DATA 

SUBROUTINES USED NAMED (CODE. KTESXi COLOR t SERCHt MATRL. TAPESi JULDYt 
WATLV) 

DEFINE FILE 1 ( 1000 1 161 .U , I SEC } 2 I 21 1 1000 U , ISECl ) 3 ( 10 .320.U IM) 
DIMENSION V(21 ) .VVI21) .VARI8) 

DIMENSION MAJTSI 10) tMACQLI 10 ) .WAT EL { 1000 ) i A INC< 20 ) . IPUNHI6) 
COMMON I WELL (1000) IVAR(8) .TVALUI81 ilADV(8l iWVARUOOQ) ISTORI2U 
COMMON HEAD(2i252) iTITLE I 2 i 108 I i JWELL t 1000 I iIADVR<8) t STRATI 208> 
COMMON ISECl ISECltISEC2A<VAL(10> iSVAL(lO) KTEST KSAT .LTEST 
HEADERS READ IN FOR LABELING VARIALSLES. 
CONTROLS FOR SELECTING WELL VARIABLES. 
FORMAT< 1615) 

NWELL-NUMBER OF WELLS TO BE SEARCHED. 

LVAR -NUMBER OF LAYER VARIABLES TO BE TESTED. 

LVAR1-NUMBER OF LAYER TEST VALUES. TVALU (IJi TO BE READ IN. 

VALUES OF (TVALU (I) FORMAT ( 6F12 .4 t /2F12 .4 ) ) 

NUM -NUMBER OF ADDITIONAL LAYER VARIABLES NOT USED FOR SEARCHING. 

NVAR -NUMBER OF WELL VARIABLES 

ICOND-'O'-BYPASS 

'1 '-CALLS IN SUBROUTINE FOR TESTING TYPE OF MATERIAL IN THE LAYERS, 
ICOL- 'O'-BY-PASS 

'I'-CALLS IN SUBROUTINE FOR TESTING LAYER COLORS. 
IDAT- 'Qt-BY-PASS 

'1 '-READS IN SUBSTRATUM ELEVATION TO BE TESTED. (SUBDA) 
INC - 'O'-BY PASS 

'N'-READS IN NUMBER OF INCREMENTS IN FEET ABOVE BEDROCK TO BE 

SEARCHED. 
I LOG- 'O'-WRITES OUT LAYER VARIABLE SEARCHED 

'I'-WRITES OUT ALL LAYERS IN A WELL 
CONTINUE 
ILVL- '0' BY PASS 

'!' WATER LEVEL CALLED FOR SEARCHING IN SATURATED ZONE 
IGSUB-'O' BY PASS 

'I 1 STEPS DOWN INCREMENTALLY BELOW GROUND SURFACE 
KTESK-'O' BY PASS 

( N' NUMBER OF PERMEABILITY VALUES TO BE USED 
KTESS-'O' BY PASS 

N' NUMBER OF STORAGE COEFFICIENT VALUES TO BE USED 
KSAT- '0' BY PASS 

, 'I 1 K AND S VALUES ARE AVERAGED FOR TOTAL THICKNESS 

PUNCH A 'I 1 IN ILVL AND A 'C 1 IN KSAT TO AVERAGE K AND S VALUES IN THE 
SATURATED ZONE 

FORMAT FOR VARIABLES NWELL THROUGH KSAT IS 1515. 
CONTINUE 

IPUNH IN) - '0' NO PUNCH '!' PUNCH 

THE FOLLOWING PUNCH COMMANDS ARE IN SIX SEPARATE FIVE COLUMN FIELDS. 
THE NUMBER FOLLOWING EACH COMMAND IS THE FIELD NUMBER. 
LAYER DATA (1) 
NET THICKNESS (2) 
WELL VARIABLES (3) 
COLOR OF LITHOLOGY 14) 
MATERIAL TYPE OF LITHOLOGY (5) 

WGT. AVE. PERM. AND WGT. AVE. STOR. COEFF. (6) 
VALUES OF SUBDA. SUSTS. SUBIN (FORMAT 13F10.0)] 
SUBDA-SUBSTRATUM ELEVATION 
SUBTS-MINIMUM ELEVATION TO BE SEARCHED 

SUBIN-INCREMENT IN FEET TO BE SUBTRACTED FROM SUBDATUM 

SUBROUTINE TAPES LOADS TEST HOLE AND LAYER DATA FROM TAPES TO DISC FILE 
FOR SEARCHING. 
CALL TAPES 
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65 


DC 


66 


C RI 


67 


1000 Ri 


68 


106 F< 


69 


II 


70 


C W 


71 


w 


72 


D 


73 


C R 


74 


999 R 


75 


107 F 


76 


I 


77 


C W 


78 


V, 


79 


1012 R 


80 


l\ 


81 


C 1 


82 


1 


83 


C F 


84 


4000 F 


85 


1 


86 


100 1 


87 


C 


88 


4001 


89 


C ! 


90 


4002 


91 


C 


92 


4003 


93 


C 


94 


4004 


95 


101 


96 


C 


97 


4005 


98 


C 


99 


4006 


100 


C 


101 


4007 


102 


C 


103 


4008 


104 


C 


105 


4009 


106 


C 


107 


3021 


108 


3024 


109 


1112 


110 




111 


3022 


112 




113 


C 


114 


C 


115 


3023 


116 


500 


117 


C 


118 


3124 


119 


C 


120 


5001 


121 


C 


122 


5002 


123 


5003 


124 


5004 


125 




126 




127 


1120 


128 


1121 


129 


1122 


130 


1123 


131 


1124 


132 




133 




134 




135 


5005 



1000 J=l2 

D LAYER VARIABLE NAMES FROM CONTROL CARDS 
READ(5106) iHEADUtll i I =1.108) 
FORMAT<72Al/i36Al) 

1 

WRITES LAYER VARIABLE NAMES ON DISC 
WRITE13 1 IMJHEAD 

999 J=l2 

READ TEST HOLE VARIABLE NAMES FROM CONTROL CARDS 
READ I 5. 10?) (HEAD(J.nIl2S2> 
FORMAT (72A1/72A1/72A1/36A1) 
1M3 

WRITES TEST HOLE VARIABLE NAMES ON DISC 
WRITEI3' IMJHEAD 

READ! 5 i 100)NWELLLVARLVAR1 tNUM .NVAR i ICOND. ICOL IDATiINCi ILOGilL 
IVL.IGSUB.KTESKiKTESS.KSAT 

TEST ON NUMBER OF TEST HOLES TO BE SEARCHED 
IF(NWELL)1013tl0134000 

READ TEST HOLE NUMBERS TO BE SEARCHED 11615) INTO ARRAY 
READ! 5 t 100) (IPUNHt I) .1=1)61 
READ(5100) I I WELL (I) t I 1 (NWELL ) 
FORMATU6I5) 

TEST ON NUMBER OF LAYER VARIABLES TO BC TESTED 
IF(LVAR)4003t40034002 

READ VARIABLES TO BE TESTED INTO ARRAY 
READ15.100) ( WARM 1 . I1LVARI 

TEST ON NUMBER OF TEST VALUES TO BE TESTED FOR 
IF (LVARl 14005.4005. 4004 

READ VALUE OF LAYER VARIABLE TO BE TESTED INTO ARRAY 
READ! 5 1 101) (TVALUU) il"liLVARU 
FORMAT l6F12.4/2Fl2i4) 

TEST ON NUMBER OF ADDITIONAL LAYER VARIABLES NEEDED! 1-7) 
IF(NUM)40Q7 I 4007.4006 

READ ADDITIONAL LAYER VARIABLES INTO ARRAY 
READlSilOO) (IADVRI I) illNUM) 

TEST ON NUMBER OF ADDITIONAL TEST HOLE VARIABLES (1-21) 
IF INVAR 14009 (4009.4008 

READ ADDITIONAL TEST HOLE VARIABLES INTO ARRAY 
READ! 5 i 100) 1 ISTORf I ) >I=1NVAR) 
TEST IF SUBDATUM ELEVATION IS TO BE SEARCHED 
IF(IDAT)3021i3021i3023 

TEST IF A HEIGHT ABOVE BEDROCK IS TO BE SEARCHED 
1F(INC)3022.3022.3024 
READ(5illl2) (AlNC(l) lltlNC) 
FORMATI16F5.1) 
GO TO 3124 
SUBDA0.0 
GO TO 3124 

READ FROM ONE CARD.UIELEV TO BE SEARCHED < 2 ) ELEV TO STOP ON i ( 3 ) 
INCREMENT TO STEP DOWN-FEET 
READ(5i500)SUBDAiSUBTS.SUBIN 
FORMAT(3F10,0) 

TEST IF MAJOR OR MINOR LITHOLOGIC 
IF (I COND 15002. 5002, 5001 

READ CODE OF LITHOLOGY WANTED INTO ARRAY 
READ 15. 100) (MAJTS(I) I1ICOND) 

TEST IF MAJOR OR MINOR COLOR OF LITHOLOGY WANTED 
1F(ICOL)50Q4|50045003 
READ (5.100) IMACOL(I) I=>liICOL) 
LIST=NUM+LVAR 
KTESTKTESK+KTESS 
IF<KTESK)1122. 1122. 1120 
READ! 5 1121) [AKVAL1 I) I-1HTESKJ 
FORMAT! 10FS.O) 



TYPE WANTED 



READ(5.1121HSVALU) ,1-l.KTESSI 

CALL CODE(NWELL.LVARiLVARliNUM.NVARiICO.NDICOL. I DAT. INC. ILOG.IPUMH 
lILVLiMAJTSiMACOLtLlSTtSUBDA.SUBTS(SUBIN.ISTORIADVR.KSAT.IGSUB.KT 
2ESK.KTESSiAKVALtSVALiAINC) 

IFUCONDI50065006.5005 

CALL MATRL ( I COND NWELL .MA JTS 1 1 PUNH < 5 I ) 
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136 


5006 


137 


5007 


138 


5080 


139 


1 


140 




141 




U2 


2 


143 


8000 


144 


8001 


145 


8002 


146 


3014 


147 


C 


US 




149 




150 




151 


C 


152 




153 


902 


154 


C 


195 




156 




157 




158 




159 


1001 


160 


1002 


161 


C 


162 




163 


903 


164 


C 


165 




166 




167 




168 




16V 


2001 


170 


2002 


171 


1003 


172 


C 


173 




174 


5016 


175 


C 


176 


C 


177 


5017 


178 




179 


5018 


180 




181 


3012 


182 


C 


ins 


C 


184 


109 


189 




186 


108 


107 


301? 


1B8 




189 


C 


190 


1004 


191 


115 


192 




193 




195 


C 


196 


C 


197 


7001 


198 




199 




200 


1110 


201 




202 


1113 


203 


1114 


204 




205 




206 





IF( J COL I 5080.5000*5007 

CALL COLOR ( I COL *NWELL MACOL , I PUNH ( M ) 

IFI IGSUBI2.2.1 

SURDA=1.0 

AMIN"9999. 

GO TO 8002 

IF ISUBOA IB 000 t 80008002 

IFt ILVLI 8002.8002*8001 

CALL WATLV(NWELL*IWELLiWATEL) 

lF(LISm004.6004i30U 

IM1 

READS HEADERS OF LAYER VARIABLES FROM DISC 

RfAOO' IMIHEAD 

DO 1003 Kli2 

J = 

TEST FOR NUMBER OF VARIABLES TO BE SELECTED 

IF<LVAR)903i903t902 

DO 1002 JC1,LVAR 

SELECTS HEADERS FOR VARIABLES 

ISTRT"! IVARl JC)*12-11) 

IEND-1STRT+11 

DO 1001 L=ISTRT( IENO 



TITLEtKf J)HEADlKiL) 

CONTINUE 

TEST FOR ADDITIONAL LAYER VARIABLES 

IF(NUM)1003, 1003, 903 

DO 2002 JCM iNUM 

SELECTS HEADERS FOR LAYER VARIABLES 

ISTRTMIADVR(JC)*12"11) 

IEND-ISTRT+11 

DO 2001 LISTRT.IENO 



TITLC(KJ)*>HEAD<K.L) 

CONTINUE 

CONTINUE 

TtST ON SUBDATUM ELEV 

IF(SUBDA)5016. 5016*3012 

IFt ILVL3013. 3013. 5017 

LABELS TEST HOLE NO. i SUBDATUM* DEPTH TO 

TOP OF A LAYER* AND ELEVATION OF THE TOP 

WRITEU.109I 



TOP OF LAYER 
OF SATURATED 



ELEVATION 
ZONE. 



OF THE 



FORMAT! + ' WN SUB LD L ELEV SATURATED ZONE') 

GO TO 3013 

WRITEI6*109) 

LABELS TEST HOLE NO.* SUBDATUM. DEPTH TO TOP OF LAYER, ELEV OF TOP OF 

LAYER IN WHICH SUBDATUM 15 FOUND 

FORMAT! 'O 1 ) 

WRITE(6108) 

FORMAT! '+'* WN SUB LD L ELEV 1 ) 

1END=J 

DO 1004 Kl*2 

WRITES HEADERS FOR VARIABLES 

WRITE(6ill5) <TITLE<K*J) tJnltlEND) 

FORMAT ( ' ' 22Xt 96A1 I 

ELV10*0 

ICONT"0 

JCONT-0 

SET BRANCH "1 

READ A TEST HOLE NUMBER TO BE SEARCHED 

SUB-SUBDA 

ICONTMCONT+1 

IF 



GO TO 1114 
IINCINC 

DO 1111 II1IINC 
DO 10 ICltNWELL 
LTEST-0 
IBRN-1 
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207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

216 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 



5000 
502 

3001 



3002 



3003 



53 



9001 



3004 
9999 

9002 

3005 



3006 
2000 



2003 
2004 

2005 

2006 

12001 



7005 
7006 

200 
202 

104 



203 
5009 
9997 



5 
5008 



ISEC=21 

JBRN2 

SUBDASUB 

IF (SUBDA) 5000,5000,3001 

1F(ILVL)300130015Q2 

JBRN=1 

READS TOTAL NUMBER OF LAYERS FOR EACH TEST HOLE FROM DISC 

READ (2'ISEOWVAR 

IKKMWELUIC) 

BRANCHES TO DIFFERENT PARTS OF PROGRAM 

GO T0<3002.3003>53) ,IBRN 

SETS MAXLA EQUAL TO TOTAL NO. OF LAYERS FOR EACH TEST HOLE 

MAXLA=WVAR(IKK) 

GO TO 3004 

READS GROUND SURFACE ELEV FOR EACH TEST HOLE FROM DISC 

GSELV*WVAR( IKK) 

GSURFaGSELV 

SETS SECTOR TO FIND BEDROCK ELEV 

ISEC-7 

SET BRANCH 3 

IBRN=3 

GO TO 3001 

READS BEDROCK ELEV FOR EACH TEST HOLE FROM DISC 

BEDRKWVAR( IKK) 

TEST IF SUBDATUM ABOVE BEDROCK Is DESIRED 

IF t I NO 3006 1 3006, 9001 

ESTABLISH SUBDATUM ELEV ABOVE BEDROCK 

SUBDA=BEDRK+AINCU1 ) 

GO TO 3006 

TEST IF SUBDATUM ABOVE BEDROCK EXISTS 

IF(SUBDA)9999, 9999, 3005 

IF I1LVL 19002 .9002 t3005 

TEST IF SUBDATUM ASOVE BEDROCK IS DESIRED 

IF(INC)3006, 3006. 3005 

SETS SECTOR TO SELECT GROUND SURFACE ELEV. 

ISEC-6 

SET BRANCH 2 

!BRN2 

GO TO 3001 

IF(IGSU8)1200l12001,2000 

SUBDA=GSURF 

SUB.DA = 5UBDA-IIGSUBMCONT) 

SUBTS*BEDRK 

IF (SUBTS } 2003, 2003 ,2005 

WRITE16.2004MWELLUC) 

FORMATt '0' i25( '*') ,'NO BEDROCK ELEVATION FOR WELL ' , 14 ,2 5 ( ' * ' 1 , / / ) 

GO TO 10 

lF(BeORK-AMlN200612001il2001 

AMIN*8EDRK 

ISEC2IWELL( 1C) 

ILAY-1 

READ ( 1 ' I SEC2 ) STRAT 

SUM=0.0 

TH1K*0,0 

IWLLMWELLI 1C) 

TEST IF TEST HOLE HAS A LOG 

1F(MAXLAHO,10,200 

TEST HOLES WITH MORE THAN TWENTY LAYERS ARE WRITTEN ON PRINTER. 

IFIMAXLA-20)203203,2Q2 

WRIT(6104)IWELL( 1C) 

FORMAT) ' 't 'LAYER NUMBER GREATER THAN 20 ON TEST HOLE 'U) 

GO TO 10 

READS THROUGH TEST HOLE LAYERS 

IF(ILVL)9997, 9997. 5009 



DO 90 JLAY 1 iMAXLA 

SUMS LAYER THICKNESS 

SUM=SUM + STRAT( ILAY2) 

TEST IF SUBDATUM EUEV IS READ 

CONTINUE 

IF(SUBDA)3008, 3008, 3007 



270 3007 IF(5UBDA-GSELVJ40l6. 4016,10 

279 4016 IFtSUBDA-BEDRK) 10,4017.4017 

280 C SUBTRACT LAYER THICKNESS OFF GROUND SURFACE 

281 4017 ELEV=GSELV-STRAT(lLAY ,2 ) 

282 C TEST NEW GROUND SURF AGAINST SUBDATUM ELEV 

283 IF ( ELEV-SUBDA)3Q10. 3010.3009 

284 C RESFTS GROUND SURFACE TO TOP OF NEXT LAYER 
265 3009 GSELV=ELEV 

286 C INCREASE LAYER NO. BY 1 

287 1LAYMLAY+1 

288 GO TO 90 

289 C COMPUTES DEPTH TO TOP OF LAYER CONTAINING SUBOATUM 

290 3010 GO TO < 5010 5012 ) JBRN 

291 5010 JBRN=2 

292 SUM=GSURF-SUBDA 

293 STRATt I LAY . 2 ) =SUBDA-ELEV 

294 ELEV=SUBDA 

295 GSELV=SUBDA 

296 SUM=SUM+STRAT< I LAY, 2) 

297 IF (KSAT) 3007.3007. 7008 

296 7008 CALL KTESX ( I LAY .57 S8 i AKV1 1 5V 1 I 

299 5012 DTOLA=GSURF-GSELV 

300 C TEST IF ADDITIONAL LAYER VARIABLES 

301 3020 1F(NUM>301630163017 

302 3017 DO 3015 1JJ=1,NUM 

303 C SETS fc TO ADDITIONAL LAYER VARIABLE 

304 KIADVR(IJJ> 

305 C SETS IV TO THE NUMBER OF THE ADDITIONAL LAYER VARIABLES PLUS THE NUMBER OF 

306 C LAYER VARIABLES TO BE SELECTED 

307 1V-IJJ+LVAR 

306 C STORES ADDITIONAL LAYER VARIABLES IN ARRAY 

309 3015 VAR( IV)=STRAT< ILAYiKI 

310 IF (LVAR 130 16. 3025. 301 6 

311 3016 DO 3011 IJ=1LVAR 

312 C SETS K TO LAYER VARIABLE TO BE SELECTED 

313 K1VAR(IJ) 

314 C STORES LAYER VARIABLES SELECTED IN ARRAY 

315 3011 VAR( IJ)=STRAT< ILAYK) 

316 C COMPUTES DEPTH TO TOP OF LAYER WHICH CONTAINS SUBDATUM 

317 3025 ELV*GSURF-(SUM-STRAT( ILAY.2H 

318 C SETS IDT TO DEPTH TO LAYER 

319 IDT=DTOLA 

320 C SETS ISUB TO SUBDATUM ELEVATION 

321 ISUB=SU8DA 

322 IF (KSAT) 300, 300 i 7003 

323 C TEST IF LOG OF ALL LAYERS IS WANTED 

324 300 1FULOG17050, 7050, 7003 

325 C WRITE DATA OUT FOR ALL LAYERS 

326 7050 IFUGSUB)7051t7051.7003 

327 C LABELS SUBDATUM ELEVATION, DEPTH TO LAYER. ELEVATION OF TOP OF LAYER 

328 C CONTAINING SUBDATUM, AND SELECTED VARIABLES. 

329 7051 IF(ELV-ELV1)7003,10.7003 

330 7003 WRITE (6.113) IWLL. ISUBi IDT i ELV (VARU I- l.L 1ST ) 

331 113 FORMATl 1 ' i I 3 I 5 14 FS. 2 t 8F 12 .4) 

332 C TEST FOR PUNCHED OUTPUT. 

333 IF(IPUNHIl) 194.94,1050 

334 1050 WRITEI7.1051) I WLL t I SUD , I DT t EL V, I VAR t I) I = 1 tLlST) 

335 1051 FORMAT! 13. I 5 I 4 tF8 .2 5F12 .4 / 5F1 2*4 ) 

336 94 1F(KSAT>301, 301.93 

337 301 IFULVLUO.10.95 

338 95 ILAYILAY+1 

339 C TEST LAYER VARIABLE VALUES SEARCHED. 

*40 C LVAR1 IS THE NUMBER OF VALUES SELECTED FOR EACH LAYER VARIABLL TESTED. 

341 3008 IFILVAR1)4018)4018.904 

342 904 00 6 JJ*=1LVAR 

343 K=IVAR(JJ) 

344 DO 6 J=1LVAR1 

345 C TEST DATA IN LAYER AGAINST TEST VALUE 

346 IFlTVALUt J) -STRATI JLAY*K1 16,4011,6 

347 6 CONTINUE 

348 GO TO 9 
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349 


C 1 


350 


4011 1 


351 


4018 1 


352 


4010 C 


353 


\ 


354 


C I 


355 




356 


8 i 


357 


C 


358 


6003 


359 


6000 1 


360 


C 


361 




362 


C 


363 


7 


364 


C 


365 


6001 


366 


6002 


367 


6020 


368 


103 


369 


C 


370 




371 


6006 


372 


6007 


373 


C 


374 


9 


375 


6008 


376 


91 


377 




378 


92 


379 


C 


380 


C 


381 


920 


382 


90 


383 


93 


384 


6009 


385 




386 




387 


6010 


388 




389 




390 


11100 


391 


11101 


392 


C 


393 


6011 


394 


4019 


395 


161 


396 




397 


C 


398 


4020 


399 


10 


400 


nil 


401 


c 


402 


6005 


403 


98 


404 


99 


405 


C 


406 


7000 


407 


C 


408 


9003 


409 


9004 


410 




411 


7002 


412 




413 


C 


414 


9005 


415 


C 


416 


C 


417 


6004 


418 


C 


419 
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TEST FOR ADDITIONAL LAYER VARIABLES 

THIK=THIK+STRAT(JLAY.2) 

IFINUM16003(6003(4Q10 

00 8 IJ=1(NUM 

K=IAOVR[ IJ) 

SUMS NUMBER OF LAYER VARIABLES TESTED AND ADDITIONAL LAYER VARIABLES 

IV=IJ+LVAR 

VAR( IVJ-STRATt JLAYiK) 

TEST ON NUMBER OF LAYER VARIABLES TO BE SEARCHED 

IF(LVAR)6001,60Qlt6000 

DO 7 IJ=1(LVAR 

SETS K TO THE VARIABLE TESTED 

KMVAR(IJ) 

SETS VARUJ) TO THE VALUE STORED IN LAYER DATA 

VAR{ IJ)*STRAT( JLAY.K) 

TEST ON VARIABLES TO BE SEARCHED FOR OR ADDITIONAL LAYER VARIABLES 

IFILIST)60046004|6002 

IFI ILVL)6020,6020,9 

WRITE! 6 (103) IWLLJLAYt(VAR{ I ) 1*1 iL 1 ST } 

FORMAT) ' ' >10X(2I5.8F12.4) 

TEST FOR PUNCHED OUTPUT 

IF(IPUNHIl) 19(9(6006 

WRITEI7i6007) IWLLtJLAYt (VAR{ IJ)iUlLIST 

FORMAT!2I4>8F9*2/(8X8F9.2/(8Xi5F9.2) 

RETURN TO READ NEXT LAYER 

lF(KTESmO(90(6008 

lFULVL)92i92i91 

KLAY=ILAY-1 

GO TO 920 

KLAY-JLAY 

SUBROUTINE KTESX COMPUTES AVERAGE PERMEABILITY AND STORAGE COEFFICIENT 

VALUES. 

CALL KTESX<KLAY(S7iS8tAKVlp5Vl) 

CONTINUE 

IFIKTEST16011(6011(6009 

AKV1=AKV1/S7 

SV1=SV1/S8 

WRITE(6(6010)IWLL,AKV1SV1 

FORMAT! 1 ( 'HOLE NUMBER= '(I4,/,' WOT. AVEPERM ( 
1AVE. STORAGE COEFFICIENT* 'iFlO.3) 

1FUPUNHI6M 6011 (6011. 11 100 

WRITE! 7tlll01 I IWLL(AKV1 (SV1 

FORMAT! 14(F5,0F5.3) 

RETURN TO READ NEXT TEST HOLE 

IF(LVARlHOtlO)4019 

WRlTEl6il61MWLLTHIK[TVALU< J) J"1LVAR1) 

FORMAT! 1 '' HOLE NO. SI5* 1 NET THICK ( ( F12.0' CODES '5X6F4.0) 

IF!IPUNH!2))10ilOi^020 

PUNCHES DATA ON CARDS 

WRITE! 7 1 161) IWU(THIK(TVALU(J) .J = 1(LVAR1) 

CONTINUE 

CONTINUE 

TEST IF SUBDATUM ELEVATION IS SELECTED 

IF(ILVL)99(99(98 

SUBDA=0.0 

I FISUBDA) 6004(600^,7000 

TEST IF POINT ABOVE BEDROCK IS SELECTED 

IFIINC)9003(9003(600< 

STEP DOWN SUBDATUM TO BE SEARCHED 

IF (IGSUB) 9004(9004.9005 

SU8DA=SUBDA~SUBIN 

IF I5UBDA-SUBTSJ6004 i 6004, 7002 

ELV1ELV 

GO TO 7001 

TEST IF SUBDATUM WANTED IS BELOW LOWER LIMIT 

IF (SUBDA-AMIN 16004, 6004,7001 

RESET ELFVATION OF TOP OF LAYER 

TEST IF ADDITIONAL TEST HOLE VARIABLES ARE SELECTED 

IF|NVAR)1011,1011(11 

SET DISC SECTOR 



l F10.0(/, t WGT, 
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420 C READ HEADERS FOR ADDITIONAL TEST HOLE VARIABLES 

421 READ <3' IM) HEAD 

422 JSTRT=1 

423 JEND=9 

424 DO 1009 IL=1*3 

425 C TEST IF NUMBER IS MORE THAN NINE 

426 IF (NVAR-JEND) 1010,1010.12 

427 C SETS NUMBER OF VARIABLES TO BE WRITTEN 

428 1010 JEND=NVAR 

429 12 DO 1007 K=l2 

430 J=0 

431 C READS NINE HEADERS FROM DISC STORAGE 

432 DO 1006 JC=JSTRTtJEND 

433 ISTRT=( I STOR( JC 1*12-11) 

434 IEND-ISTRT+11 

435 DO 1005 L=ISTRTIEND 

436 J*J+1 

437 1005 TITLE<KJ)=HEADIKL) 

438 1006 CONTINUE 

439 1007 CONTINUE 

440 IEND=J 

441 DO 1008 K=lt2 

442 C WRITES NINE HEADERS ON EACH PRINTER LINE 

443 1008 WRITE (6 ( 105) (TITLE I K . J ) J= 1 1 1 END ) 

444 105 FORMAT) ' ' ilOSAl) 

445 C TEST IF SUBDATUM ELEVATION IS SELECTED 

446 DO 22 J=1NWELL 

447 00 20 I=JSTRTJEND 

448 C STORES SELECTED VARIABLES ON DISC 

449 ISEC1=ISTOR! I) 

450 C READS VARIABLES FROM DISC 

451 READ (2'ISECDWVAR 

452 C SETS IK TO TEST HOLE NUMBER 

453 IK*IWELL(J) 

454 C SETS IKK TO VARIABLE 

455 IKKMSTOR(I) 

456 V< IKK)=WVAR< IK) 

457 20 CONTINUE 

458 DO 21 I=JSTRT.JEND 

459 KMSTOR(l) 

460 21 VV( I )=V(K) 

461 C WRITES VARIABLES FOUND AT SUBDATUM ON PRINTER 

462 WRITE (6i5400) ( VV ( I) i I = JSTRT t JEND ) 

463 5400 FORMATf' ( i9F12.3) 

464 C TEST FOR PUNCHED OUTPUT 

465 IFUPUNHO) )22221054 

466 1054 WRITEI7.1055MVVU) I "J5TRT i JGND ) 

467 1055 FORMAT < 6F12* 3 / 6F12 . 3 t / i 6F12 .3 . / 3F 12. 3 ) 

468 C PUNCHES VARIABLES FOUND AT SUBDATUMi SIX PER CARD 

469 22 CONTINUE 

470 C TEST IF ALL VARIABLES SELECTED HAVE BEEN WRITTEN 

471 IF ( JEND-NVAR)23tl0111011 

472 C RESET TO NEXT NINE HEADINGS AND VARIABLES 

473 23 JSTRT=JEND+1 

474 JEND=JEND+9 

475 1009 CONTINUE 

476 1011 CONTINUE 

477 GO TO 1012 

478 1013 STOP 

479 END 

480 C 

481 C 

482 C SUBROUTINE CODE DISPLAYS A LIST OF OPTIONS SELECTED BY THE USER. 

483 C RETRIEVAL CONTROL CARD CODING 15 DISPLAYED ON PRINTER FOR EACH JOB SETUP. 

484 SUBROUTINE CODE ( NWELL iLVAR tLVARl NUM tNVAR . 1 COND 1 1 COL 1 1 OAT t INCt I LOG 
4ft 5 llPUNHMLVLiMAJTSMACOLlLlSTiSU8DAfSUBTStSUBINtlSTORIADVKtKSATtl 

486 2GSUB.KTESKKTESS.AKVALSVAL.AINC) 

487 DIMENSION AKVALf 10 ) tSVAL ( 10 ) .A INC (20) 

488 DIMENSION MAJTS! 10) tMACOLtlO) .ISTORI21) iIADVR(8J IPUNH16) 

489 COMMON IWELL I 10001 t IVARI 8 I TVAUI( 8 I 

490 WRITEI6.120) 
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492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 



120 FORMAT) 'I 1 i 1 IF TEST CODES ARE EQUAL TO 1, THEN TEST IS USED. IF E 
10UAL TO Oi THEN TEST IS BY-PASSED '.// 1 
WRITE!6100) NWELL 

100 FORMAT! 1 '. 'NUMBER OF TEST HOLES TO RE SFARCHED= 'iI3l 
WRITE(6 ( 133) { I WELL! I I il=l .NWELL) 

133 FORMAT! 1 '('TEST HOLES THAT ARE BEING TESTED ARE ' . / .41 I 24 I 5 . / I ) 
WRITE(6.10ULVAR 

101 FORMAT! '0' t 'NUMBER OF LAYER VARIABLES TO BE TE5TflD = ',12) 
IF ILVAR)22.1 

WRITE! 6il02) { WARt I) iIMtLVAR) 



1 
102 

2 
103 



FORMAT! 1 '('THE VARIABLES THAT ARE TO BE TESTED ARE'/, 815) 

WRITEI6.1031LVAR1 

FORMAT! '0' ' NUMBER OF TEST VALUES* '(12) 



IFILVAR1)443 

3 WRITEI6il04) (TVALUI I) I=1.LVAR1I 

104 FORMAT!' 't'TEST VALUES OF EACH LAYER ARE ' > / t 8F10*4 } 

4 WRITF(6(105INUM 

105 FORMAT! '0' t 'NUMBER OF ADDITIONAL LAYER VARIABLES= 't!2) 
IF(NUM)6.6t5 

5 WRITE!6d06) 1 IADVR! I I t I=1NUM) 

106 FORMAT! ' '('ADDITIONAL LAYER VARIA8LES ARE'/(8I5) 

6 WRITE!6,116)ILIST 

116 FORMAT! '0' t'TOTAL NUMBER OF LAYER VARIABLES THAT ARE CALLED= P I2) 
WR!TE!6il07)NVAR 

107 FORMAT! '0' I'NUMBER OF TEST HOLE VARIABLES* '(12) 



7 
108 



IF!NVAR)8i8t7 

WR1TEI6108) ( ISTORl I) . I=1.NVAR) 



FORMAT! 



.'THE TEST HOLE VARIABLES ARE't/21I5) 



8 WRITE(6il09HOAT 

109 FORMAT! '0' . 'SUBDATUM ELEVATION CODE= SI2) 
IF!IDAT)10109 

9 WRITE<6110)SUBDA.5UBT5.SUBIN 

110 FORMAT! 1 '.'ELEVATION TO BC SEARCHED* '.F10.0(/,' ELEVATION TO TER 
1MINATE SEARCHING" 'F10.0/' STEP-DOWN INCREMENT*: SF10.0) 

10 IF!INC)l32il32.130 

130 WRITE16.131MNC. IAINC! I) ( 1 = 1 (INC I 

131 FORMAT! 'ONUMBER OF HEIGHT INCREMENTS ABOVE BEDROCK TO BE SEARCHED 
1= '.I2./,' THE INCREMENT VALUES ARE ' i / t 20F1 0. 2 I 

132 WRITE(6ill2)lCOND 

112 FORMAT! 'O 1 . 'TEST IF MAJOR OR MINOR LITHOLOGY TYPE WANTED* ',121 
IF(ICOND)12. 12.11 
WRITE(6dl3l (MAJTS! I ) . I=1.ICONDI 



.'LITHOLOGY TYPE CODES ARE'i/tBIB) 



IF MAJOR OR MINOR COLOR OF LITHOLOGY WANTED" 1 I3 



11 

113 FORMAT! 

12 WRITE(6.114) ICOL 

114 FORMAT! '0' t 'TEST 
IF! ICOL)1414)13 

WRITE! 6. 115) (MACOL! I) .1=1 ( ICOL) 

FORMAT!' 't'LITHOLOGY COLOR CODES ARE',/, 815) 

WRITE!6,117)ILVL 

FORMAT! '0' ( 'TEST IF WATERLtVEL OF SATURATED ZONE IS TO BE USED* ' 
112) 

WRITE!6(12l)KSAT 
121 FORMAT! 'OTEST IF SATURATED ZONE IS TO BE LOGGEP= P ,I2) 

WRITE(6lie)ILOG 

FORMAT! '0' t'TEST IF ALL LAYERS ARE TO BE LOGGED= '.12) 

IF(IGSUB)124.124,122 

WRITE(6tl23)IGSUB 

FORMAT! 'OSTEP DOWN INCREMENT FROM GROUND SURFACE ELEVATION* '.I2( 
1* FEET 1 I 

124 IF(KTESK)129,129.125 

125 WRITE!6il26)KTESK 

126 FORMAT! 'ONUMBER OF K-PERMEABI LI TY AND S-STORAGE COEF. VALUES* ',12 
1) 

WRITE!6,127MAKVAL(I) 

127 FORMAT! 



13 
115 

14 
117 



118 

122 

123 



i I 



('THE <-PERMEABILITY VALUES ARE ' /i 10F8.0) 



WRITE(6(128MSVALII)*I1KTSS) 

128 FORMAT!' ' . ' THE S-STORAGE COEF. VALUES ARE ' / . 10F8 .3 I 

129 WR!TE!6.119) (IPUNH! I I tl = l .61 

119 FORMAT! '0' ( 'TEST FOR PUNCHING DATA '// 5X (' LAYER DATa= l tI2, 
1NET THICKNESS^ ' 1 1 2 ( / t 5X ' TEST HOLE OATA= ' . I 2 ( / 5X ' COLOR OF 
1 ' (I2(/(5X( 'MATERIAL TYPE OF LITHOLOGY* SlSi/tSXi 1 WGT. AVE 



LITHOLOGY" 
PFRM 
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562 2. AND WGT. AVE. STOR. COEF.= '.12) 

563 RETURN 

564 END 

565 C 

566 C 

567 C SUBROUTINE KTESX COMPUTES AVERAGE K. AND S VALUES WEIGHTED ON LAYER 
<J68 C THICKNESS 

569 SUBROUTINE KTESX ( I LAY S7 58 AKV1 iSVl ) 

570 COMMON I WELL (1000 I IVAR (8) tTVALU(B) . IADV18I WVAR(1000) iISTORtM I 

571 COMMON HE AD(2. 252). TITLE(2108)iJ WELL (1000) IADVR(8) i STRATI 20*6) 

572 COMMON I SEC , I SECl I SEC2 tAKVAL(lO) .SVALl 10 I KTEST KSAT . LTEST 

573 IFILTESTJ10. 10*15 

574 10 S7*0.0 

575 S8-0.0 

576 AKV1=0.0 

577 SVlnQ.O 

578 LTEST=LTEST+1 

579 15 J7=STRAT( ILAY.7) 

580 J8=STRAT( ILAY.6) 

581 IF(J7)2121.20 

582 20 S7=S7+STRAT(ILAY.2) 

583 AKV1 = AKV1-HAKVAL<J7)*STRAT( ILAY.2H 
504 21 IFIJ8126. 26,25 

585 25 S8=S8+STRAT< ILAYiZI 

586 SV1=SV1+(SVAL( J8)*STRATI ILAYtZ) I 

587 26 RETURN 

588 END 

589 C 

590 C 

591 SUBROUTINE COLOR ( ICOLR .NWELL .MACOLt 1 PUNH ) 

592 C SUBROUTINE TO READ AND TEST FOR MAJOR OR MINOR COLOR OF LITHOLOGY 

593 DIMENSION MACOL(IO) 

594 COMMON I WELL (1000) . WAR ( 8 1 iTVALUlfll ( IADVI8] tWVARUOOO) i!STOR(21) 

595 COMMON HEAD ( 2 2 52 ) i T I TLE ( 2 1 108 ) JWELL ( 1000) iIADVROJ tSTRAT1208) 

596 COMMON I SEC i I SEC 1 1 1 SEC2 

597 100 FORMAT! 1 '.'HOLE NO. ' t 2X t* LAYER NO. ' 2X i ' MAJOR COLOR ' . 2X i ' M I NOR CO 

598 1LOR' .4X'ELEV ' i 3X ,' DEPTH ' i 3X TH ICKNESS ' i 5X ( NET TH ICK < MAJOR I NE 

599 2T THICKIMINOR) " J 

600 101 FORMAT!' ' 2X i I 3 i 6X I 2 i 10X i 12 20X tF6 1 i 2XF6 1 .6X F6. 2 ) 

601 102 FORMAT!' ' 2X i I 3 8X I 2 *23X i I 2 i 7X i F6 . 1 i 2X iF6. 1 6X F6. 2 > 

602 103 FORMAT! '+' .BOX.F6.2) 

603 104 FORMAT! '+' .99X.F6. 2) 

604 C WRITE HEADERS FOR COLORS FOUND IN LlTHOLOGY 

605 WRITE(6.100) 

606 C ICOLR IS THE NUMBER OF COLORS SELECTED 

607 00 6 K=1ICOLR 

608 C MACOL(K) IS THE NUMERICAL CODE FOR EACH COLOR SELECTED 

609 JTEST=MACOL(K) 

610 C NWELL IS THE NUMBER OF TEST HOLES TO BE SEARCHED 

611 00 5 IC1*NWELL 

612 C CALLS SUBROUTINE NAMED 'SERCH 1 

613 CALL SRCHINLAYGSELVIC) 

614 C SETS SECTOR TO TCST HOLE NO. BEING SEARCHED 

615 ISEC2=IWELL(IC) 

616 C READS FROM DISC LAYER DATA 

617 READ (l'ISEC2ISTRAT 

618 SUM=0.0 

619 SUM2=0.0 

620 SUM3=0.0 

621 J=l 

622 C SETS IWLL TO TEST HOLE NO. BEING SEARCHED 

623 52 IWLL=IWELL(IC) 

624 C READS LAYER DATA FOR ALL LAYERS IN EACH TEST HOLE 

625 53 DO 4 JJMiNLAY 

626 C SETS TYPE TO COLOR OF LAYER 

627 50 TYPE=STRAT!J.4> 

628 C SETS THICK TO THICKNESS OF LAYER BEIMG SEARCHED 

629 THICKSTRAT(J.2) 

630 C SUMS LAYER THICKNESS OF ALL LAYERS AS SEARCHED 

631 SUM=SUM+STRAT(Jt2) 

632 C KTEST IS A THREE DIGIT NUMBER [COLOR CODE) OF LAYER 
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633 KTEST= IF IX ( TYPE ) /1Q 

634 C TEST IF COLOR IN LAYER IS SELECTED COLOR 

635 IF (KTEST-JTEST) 2. 1. 2 

636 C SDMI IS THE THICKNESS TO TOP OF LAYER WHERE COLOR IS FOUND 

637 1 SUM1 = SUM-STRATU2) 
636 5UM2-SUM2+TH1CK 

639 C COMPUTES ELEVATION OF TOP OF LAYER WHERE COLOR IS FOUND 

640 ELEV=GSELV-SUM1 

641 C ICOL IS TWO NUMBER CODE 

642 ICOL TYPE 

643 C LABELS HOLE NO. i LAYER NO.* COLOR CODE, ELEV OF TOP OF LAYERt DEPTH TO 

644 C LAYERt LAYER THICKNESS 

645 WR1TE(6.101IIWLL J J . I COL .ELEV SUM1 i TH I CK 

646 C TEST FOR PUNCHED OUTPUT 

647 IFUPUNHI4i4|10 

64B C PUNCHES HOLE NO.. LAYER NO.. COLOR CODE ELEV. OF TOP OF LAYER. DEPTH 

649 C TO LAYER* LAYER THICKNESS ON ONE CARD 

650 10 WRITE<7,11HWLL,JJICOLELEVSUM1, THICK 

651 11 FORMAT) 1 COLORl ' .3 15 F6 . 1 .2X tF6. 1 . 6X.F6. 2 

652 GO TO 4 

653 2 KTEST=TYPE-FLOAT1KTEST)*10. 

6**4 C TEST IF COLOR ASKED FOR IS MINOR COLOR 

655 IF(KTEST-JTESTI4i3i4 

656 C CALCULATES DEPTH TO TOP OF LAYER 

657 3 SUM1SUM-STRAT(J2) 

658 SUM3SUM3+THICK 

659 C CALCULATES ELEV. TOP OF LAYER 

660 ELV=GSELV-SUM1 

661 C ICOL IS A TWO NUMBER CODE 

662 ICOL'TYPE 

663 C LABELS HOLE NO.. LAYER NO.i COLOR CODE, ELEV OF TOP OF LAYER, DEPTH 

664 C TO LAYER. LAYER THICKNESS 

665 WRITE (6. 102) IWLL, JJ t ICOL .ELEV tSUMl . TH ICK 

666 IF(IPUNH144,12 

667 C PUNCHES HOLE NO.. LAYER NO., COLOR CODE, ELEV <JF TOP OF LAYER, DEPTH 

668 C TO LAYERt LAYER THICKNESS ON ONE CARD 

669 12 WRITE(713HWLL,JJ, ICOL lELEVi SUM1 THICK 

670 13 FORMAT! 'COLOR2 ' 3 I 5 iF6 . 1 2X ,F6. 1 6X . F6. 2 ) 

671 4 J-J+1 

672 IFlSUM2)5555,5<t 

673 54 WRITE<6tl03)SUM2 

674 55 IF|SUM3)5t5,56 

675 56 WRITE16.104ISUM3 

676 5 CONTINUE 

677 6 CONTINUE 

678 RETURN 

679 END 

680 C 
661 C 

682 SUBROUTINE SERCH (NLAY tGSELV. I C ) 

683 C SUBROUTINE READS FROM DISC (II TOTAL NO. OF LAYERS. (2) GROUND 

684 C SURFACE ELEVATION FOR SELECTED TEST HOLE 

685 COMMON I WELL (1000 1 ( WAR 1 8 ) iTVALU ( 8 I ,IADVI81 iWVARtlOOO) iISTORlZl) 

686 COMMON HEAD (2, 2521 , T I TLE I 2 108 I . JWELL ( 1000) tlADVRU) , STRATI 20, 3 I 

687 COMMON ISEC , ISEC1 I SEC2 

688 C SET SECTOR NUMBER TO READ TOTAL NO. OF LAYERS FROM DISC 

689 IBRN0 

690 ISEC*21 

691 3001 READ12 1 ISECJWVAR 

692 IKKKIWELL(IC) 

693 IFUBRNI10, 10,3002 

694 10 NLAYWVAR(IKKK) 

695 C SET SECTOR NUMBER TO READ GROUND SURFACE ELEV. FROM DISC 

696 ISEC*6 

697 IBRN'l 

698 GO TO 3001 

699 C SETS GSELV TO GROUND SURFACE ELEVATION 

700 3002 GSELV=WVAR< IKKK) 

701 RETURN 

702 END 

703 C 

38 



704 


C 




705 






706 


c 




707 


C 




708 






709 






710 






711 






712 




100 


713 






714 






715 




101 


716 




102 


717 




103 


718 




104 


719 






720 


c 




721 






722 


c 




723 






724 


c 




725 






726 


c 




727 






728 






729 


c 




730 






73L 






732 






733 






734 






735 




52 


736 


c 




737 




53 


738 


c 




739 




50 


740 


c 




74L 






742 


c 




743 






744 


c 




745 






746 


c 




747 






748 






749 




30 


750 






751 






752 




31 


753 




1000 


754 




32 


755 






756 




1001 


757 




33 


758 






759 




34 


760 






761 






762 




35 


763 






764 






765 




1002 


766 


c 




767 




1 


76S 


c 




769 






770 


c 




771 






772 


c 




773 






774 


c 





SUBROUTINE MATRL < ICOND tNWELL .MAJTS i IPUNH ) 

SUBROUTINE SELECTS MAJOR AND MINOR MATERIAL TYPE OF 

LITHOLOGY 

DIMENSION MAJTS (10) 

COMMON I WELL! 1000) .IVARU) .TVALUI8) tlADVIB) tWVAR 11000) t ISTOR121) 

COMMON HEAD! 2.252) .TITLE I 2.108) . JWELLI 1000) .IADVRIS) .STRAT(20t8) 

COMMON ISEC.ISEC1.ISEC2 

FORMAT!' ' 'HOLE NO. ' . 2X i ' LA YER NO ' 2X t ' MAJOR TYPE'3X, 'MINOR TY 
1PE' .5X 'ELEV .4Xt 'DEPTH' .3X, 'THICKNESS' 5Xi 'NET THICK (MAJOR) NET 
2THICKIMINOR) ' ) 

FORMAT! ' l 2X.I38XI210XI5.18XF>.l2XF6.1i6XtF&.2) 

FORMAT! ' l 2xl3aXI222X.I5.6X(F6.1i2X.F6.1|fcXiF&.2) 

FORMAT! '+' i80XF6.2 J 

FORMATt '+ 1 99XiF6.2) 

WRITE16.100) 

ICOND IS THE NUMBER OF MATERIAL TYPES SELECTED. 

DO 6 K-l, ICOND 

SETS JTEST EQUAL TO MATERIAL CODE SELECTED <2 DIGIT) 

JTEST=MAJTS(K) 

SEARCH THROUGH EACH SELECTED TEST HOLE 

DO 5 IC=1NWELL 

CALLS SUBROUTINE SERCH 

CALL SERCH ( NLAY .GSELV. I C ) 

ISeC2=IWELL( 1C) 

READ LAYER DATA FROM DISC FOR EACH HOLE 

READ !1 ' ISEC2 JSTRAT 

SUM=00 

J=l 

SUM2=0.0 

SUM3=0.0 

IWLL=IWELL( 1C) 

SEARCH ALL LAYERS IN EACH HOLE 

DO 4 JJ=1NLAY 

SET TYPE EQUAL TO MATERIAL CODE FOR EACH LAYER (5 DIGIT! 

TYPE=STRAT( Ji3l 

SUM IS THE TOTAL THICKNESS INCLUDING THE LAYER BEING SEARCHED 

SUM=SUM+STRAT< J2) 

SETS KTEST TO A TWO DIGIT CODE OF MAJOR MATERIAL TYPE 

KTE5T=IFIXITYPE)/1000 

THICK IS THE THICKNESS OF LAYER BEING SEARCHED 

THICK=STRAT(J.2 I 

TEST IF MAJOR MATERIAL TYPE IS IN SELECTED LAYER 

IDl=TYPE/10000i 

IF(ID1-2)1002301002 

ITYPE=TYPE/1000. 

ID2=1TYPE-(ITYPE/10*10) 



ITEST=ID2+10 

IF(ITEST-JTEST)32ili32 

ITYPE=TYPE/10 

ID4=ITYPE-MTYPE/10nO) 

IF( 104-2)35.33.35 

ITEST=ID4+10 

IF( ITEST~JTEST)414 

ITYPE=TYPE/100. 

ID23=ITYPE-(]TYPE/100*100) 

IFIID23-JTESTU001, 1.1001 

ITYPE=TYPE 

!TEST=ITYPO( I TYPE/ 100*1 00) 

IF( ITEST-JTEST)4tl,4 

IF(KTEST-JTEST)2.1.2 

SUM1 IS THE DEPTH TO TOP OF LAYER WHERE MAJOR TYPE IS FOUND 

SUM1=SUM-STRAT!J.2) 

SUM2 IS THE TOTAL THICKNESS OF LAYERS WITH MAJOR TYPE 

SUM2=SUM2+THICK 

ELEV IS ELEVATION OF TOP OF LAYER CONTAINING MAJOR TYPE 

ELEV=GSELV-SUM1 

SETS TYPE TO AN INTEGER NUMBER 

ITYPE=TYPE 

LABELS HOLE NO.* LAYER NO. ( MAJOR MATERIAL TYPE. ELEV. OF TOP OF 



39 



775 


C 


LAYERt DEPTH TO LAYERt LAYER THICKNESS 




776 




WRITEI&.lOl) IWLLt J J . I TYPE. ELEV iSUMl . THI CK 




777 


C 


TEST FOR PUNCHED OUTPUT 




778 




IFI IPUNHJ44i7 




779 


C 


PUNCHES HOLE NO . LAYER N0 MAJOR MATERIAL TYPE. ELEV* 


OF TOP OF 


780 


C 


LAYER. DEPTH TO LAYER. LAY THICK 




781 




7 WRITE(7t8) IWLLi JJi I TYPE * ELEV t SUM1 (THICK 




782 




8 FORMAT! I MAJO I ?"|3162F8 < 1F10.2) 




783 




GO TO 4 




784 


C 


BREAK DOWN TYPE OF MATERIAL CODE TO TWO DIGITS 




785 




2 INUM=TYPE/100G. 




766 




ITYPE=TYPE 




787 


C 


BREAK DOWN TYPE OF MATERIAL CODE TO MINOR TYPE ONLY 




788 


C 


TEST IF MINOR TYPE OF MATERIAL CODE CALLED IS IN LAYER 




789 




INUM=TYPE-(FLOAT(INUM)*1QOO. ) 




790 




IAOJ-INUM/100 




791 




!F( IADJ-L I21i20i21 




792 




20 KTEST=INUM/10 




793 




GO TO 22 




794 




21 KTEST*INUM-(IADJ*100 




795 




22 IF(KTEST-JTEST)4.34 




796 


C 


SUML IS THE DEPTH TO TOP OF LAYER WHERE MINOR TYPE CODE 


IS FOUND 


797 




3 SUML*SUM-STRAT< J.2) 




798 




SUM3-SUM3+THICK 




799 


C 


ELEV IS THE ELEVATION OF TOP OF LAYER WHERE MINOR TYPE 


IS FOUND 


800 




ELEVGSELV-SUMl 




801 


C 


LABELS WELL NO. LAYER NO. i MINOR MATERIAL TYPE* ELEV OF TOP OF 


802 


C 


LAYER* LAYER THICK 




803 




WRITE(6il02) IWLL* JJ . I TYPE t ELEV .SUM1 THI CK 




804 




IFI IPUNHJ4.4.9 




805 




9 WRITEI710I IWLLi JJ. I TYPE . ELEV t SUM 1 .THI CK 




806 




10 FORMAT! 'M I NOR ' 3 16. 2F8. 1 F10 2 ) 




807 




4 J*J+1 




80S 


C 


TEST IF NET THICKNESS OF LAYERS WITH MAJOR MATERIAL IS 


SELECTED 


809 




IF(SUM2)5555i54 




810 


C 


WRITE NET THICKNESS OF LAYERS WITH MAJOR MATERIAL TYPE 




811 




54 WRITE(6.103)SUM2 




812 




55 IF(SUM3)5.556 




813 




56 WRITEI6.104ISUM3 




814 




5 CONTINUE 




815 




6 CONTINUE 




816 




RETURN 




817 




END 




818 


C 






819 


C 






820 




SUBROUTINE TAPES 




821 


C 


SUBROUTINE TAPES READS TAPES USED IN STORAGE AND RETRIEVAL PROGRAM 


822 


C 


TAPE NUMBER 8 FOR WELL VARIABLE DATA* 




823 


C 


TAPE DRIVE NUMBER 9 FOR LAYER DATA. 




824 




COMMON I WELL! 1000) ( I VAR ( 8 ) TVALU18) JADV(8) .WVARI 1000) . 


ISTOR(21) 


825 


COMMON HEAD{2i252liTITLE1210eiiJWELL( 1000) ilADVRf B) iSTRAT(20*8) 


826 




COMMON ISECiISECl.ISEC2 




827 




ISECl=l 




828 




DO 10 I =1.21 




829 


C 


READS TEST HOLE DATA VARIABLES FROM TAPE 




830 




READ(8END11)WVAR 




831 


C 


WRITES TEST HOLE DATA ON DISC FILE 




832 




10 WRITE12 1 ISEC11WVAR 




833 




11 DO 20 J=L.1000 




83ft 


C 


READS LAYER DATA VARIABLES FROM TAPE 




835 




READ(9iEND=30) IMiI.STRAT 




836 


C 


SELECTS LAYER VARIABLES FOR SPECIFIED TEST HOLES 




837 




1F( rM>30>3020 




8^8 


C 


WRITEUAYER DATA ON DISC) 




839 




20 WRITEilMMJSTRAT 




840 




30 RETURN 




841 




END 




842 


C 






843 


C 






844 


C 


SUBROUTINE JULDY COMPUTES JULIAN DATE FROM CALENDAR DATE AND MILITARY 


845 


C 


TIME. 





40 



846 


SUBROUTINE JULDY < IMO . I DY . I YR . IHR iT I ME ) 




847 


NYR=]YR/4 




846 


NYRNYR4 




849 


IF (NYR-IYR) 2.1*2 




850 


1 IB*1 




851 


GO TO 3 




852 


2 IB2 




853 


3 JDAY=0 




354 


00 10 Ili IMO 




855 


GOTO (4. 5.6. 5, 9. 5. 9. 5.5.9.5. 91 I 




856 


4 M0*0 




857 


GO TO 10 




858 


5 M031 




859 


GO TO 10 




860 


6 GO TO (7.8) .IB 




861 


7 M0*29 




862 


GO TO 10 




863 


8 M0=28 




864 


GO TO 10 




865 


9 M0 = 30 




866 


10 JDAY=JDAY+MO 




867 


JDAY*JDAY+IDY 




868 


DAY*JDAY-1 




869 


DAY*DAY*1440. 




870 


HR*< IHR/10Q)*60 




871 


HRlIHR-( ( IHR/100)*100> 




872 


TIME-DAY+HR+HR1 




873 


102 RETURN 




874 


END 




875 


C 




876 


C 




877 


C SUBROUTINE WATLV COMPUTES GROUND WATER ELEVATION AT SELECTED TEST HOLES 


878 


C WHEN DATA FROM THE SATURATED ZONE IS RETRIEVED 




879 


SUBROUTINE WATLV (NwELLS I WELL iWATEL ) 




B80 


C WATER LEVEL PROGRAM COMPUTES AVERAGE GROUND WATER ELEVATION FOK 


A SELECTED 


881 


C TIME PERIOD. 




882 


DIMENSION I WELL! 1000) .WATEL(IOOO) 




883 


REWIND 4 




884 


JWELL1-0 




885 


C BEGINNING (IMOi IDYt lYRi ITIM) AND ENDING (JMOi JDY JYR. JTIM) 


MONTH. 


886 


C DAY. YEAR. MILITARY TIME OF PERIOD FOR WHICH AVERAGE WATER LEVEL 


IS 


887 


C COMPUTED IS READ FROM DATA CARD 




888 


1 READ (5. 101 1 IMO.IDY. IYR . I T IM* JMO* JDY UYR t JTIM 




889 


101 FORMATU6I5) 




890 


WRIT(6110MMO.IDY.IYR,ITIM.JMQ,JDY.JYR|JTIM 




891 


110 FORMAT! '0' ,25( '*' I , "TIME PERIOD FOR SEASONAL ' i 12 i '/ ' 12 ' / ' I2i 


'/ 


892 


l'.I4, ( TO l iI2i'/ l >I2> l / l iI2i l / l iI^ilXi25P* l )i//) 




893 


1IYR-IYR 




894 


IFt IMOJ107.107.5 




895 


C SUBROUTINE JULDY COMPUTES TIME IN MINUTES FROM MONTH. DAY. YEAR. 


ANO 


896 


C MILITARY TIME 




897 


C BTIM-BEGINNING TIME (JULIAN DATE IN MINUTES) OF PERIOD FOR WHICH 


WATER 


898 


C LEVELS ARE AVERAGED 




899 


5 CALL JULDYtlMO.IDY.IYRpITIM.BTIM) 




9HQ 


JF{ IYR-JYR)2i33 




901 


C EYEAR-WHEN PERIOD FOR WHICH WATER LEVELS ARE AVERAGED CROSSES END OF YEAR 


902 


C THEN EYEAR (END OF FIRST YEAR. JULIAN DATE IN MINUTES) IS 


CALUCLATED 


903 


2 CALL JULDYtl2.31MYR.2359.EYEAR) 




904 


500 IIYR-IIYR+l 




905 


IFI I IYR-JYRJ501. 502*502 




906 


C STIME-WHEN PERIOD FOR WHICH WATER LEVELS ARE AVERAGED SPANS A FULL YEAR 


907 


C THAT YEARS TIME (JULIAN DATE IN MINUTES) IS CALCULATED 




908 


501 CALL JULDY(12.31.IIYR.2359.STIME) 




909 


EYEAR=EYEAR+STIME 




910 


GO TO 500 




911 


C EETIM-THE JULIAN DATE IN MINUTES FROM JANURARY 1 OF THE FINAL YEAR TO THE 


912 


C END OF THE SELECTED TIME PERIOD 




913 


502 CALL JULDYIJf/O.JDYJYR.JTIM.EETIM) 




914 


ETIME=EETIM 




915 


ETIM=EETIM+EYEAR 




916 


GO TO 6 





41 



917 C ET1M-END OF SELECTED PERIOD (JULIAN DATE IN MINUTES) 

918 3 CALL JULDYUMOJDYiJYRiJTlMETIM) 

919 ETIME=ETIM 

920 C TTIM-TOTAL TIME (MINUTES) OF SELECTED TIME PERIOD 
931 6 TTIMaETIM-BTIM 

922 C RATE OF WATER LEVEL CHANGE IS COMPUTED FROM ELEVATIONS READ FROM TAPE 

9?3 C WITHIN THE SELECTED TIME PERIOD. 

924 DO 1110 I=liNWELLS 

925 IFUWFLLfl ) -JWELL1 1 7 t 1000 t8 

926 7 REWIND 4 

927 JWELL1=0 

928 8 KKK=I'<JELL( I) 

929 JBRN-1 

930 JBN=1 

931 STIM=0.0 

932 SGH=0.0 

933 TIMZ=0.0 

934 TM-0.0 

935 10 READI4END107) NWELL iNCODE .NYR ,NMO iNOY iNT IMi ELEV 

936 IF (NWELL-IWFLLI I) ) 10.20.85 

937 20 CALL JULDY ( NMOiMDYiNYR MT IM.T IME ) 
939 GO T0,(3000i3001) iJBRN 

939 3000 JBRN=2 

940 IF(JYR-NYR)86t30.3010 

941 3001 IFUYR-NYRJ40. 30.3010 

942 30 IF(BT!M-TIME)3003. 3002.3010 

943 3002 WGH-0.0 

944 ELEVZ'ELCV 

945 TIMEB=TIME 

946 GO TO 3011 

947 3003 IF(NYR1~NYR)3005>30043004 

948 3004 JBN=2 

949 GO TO 40 

950 3005 JBM=3 

951 GO TO 40 

952 3010 NYR1=NYR 

953 GO TO 10 

954 40 BACKSPACE 4 

955 BACKSPACE 4 

956 READ(4tEND=107)NWELLNCODEtNYRtNMONDYtNTIMELEVl 

957 CALL JULDY 1NMO,NDYNYRMT1MT IMA 

958 45 READ(4tFND'=107)KWFLLiKCOD>KYRKMO.KDYKTIMiELEV2 

959 IFiKWELL-IWELLU ) J86i46t87 

960 46 IF((KYR-NYR)-1)3014. 3013*3012 

961 C A BREAK IN RECORD GREATER THAN ONE YEAR IS PRINTED IN OUTPUT 

962 3012 WR I TCI 6. 109) I WELL (I ] .NMO tNDY iNYR f NT iMiKMOtKDY iKYR KT IM 

963 109 FORMAT (' WELL NUMBER f il4(' BREAK IN RECORD FROM '.313*15.' TO ' 

964 1313.15) 

965 GO TO 1005 

966 3013 IF(KMO-K12-NMO)-12)3014.3012.3012 

967 3014 CALL JULDY < KMOiKDY ,KYR *KT IKi T IMEB } 

968 GO TOI206. 207. 3006) ,JBN 

969 3006 CALL JULDY ( 12 . 31 .NYR ,2359 . TM ) 

970 DIFF={ (TM-TIMA)+BTIM|-KTIMEB-STIM) 

971 DIFF1=TIMEB-BTIM 

972 GO TO 208 

973 206 CALL JULDY ( 12 .31 .NYR .2359 t TM ) 

974 DIFF=(TM-TIMA)+TIMEB 

975 DIFF1=(TM-BTIM)+TIMEB 

976 GO TO 208 

977 207 DIFF=TIMEB-TIMA 

978 DIFF1=TIMEB-BTIM 

979 208 SLOPE=(FLCV2-ELEVll/DIFF 

980 GHA*ELCV2-<SLOPE*DIFF1) 

981 DlFFaDIFFl 

982 STIM=STIMtDlFF 

983 IFISTIM-TTIM)10021011.1010 
964 1010 OIFF = DIFF-(TM-ETIM-("TIMEB) 

915 GHB=ELEV2-SLOPE*(TM-ETIM+TIMEB) 

986 GO TO 1012 

987 1011 -GH 

42 



988 


1012 


989 




990 


1002 


991 




992 


3011 


993 




994 




995 


50 


996 




997 


51 


998 




999 


200 


1000 




1001 




1002 


201 


1003 


203 


1004 




1005 


55 


1006 




1007 




lOOfl 




1009 




1010 


800 


1011 




1012 




1013 




1014 


70 


1015 


80 


1016 




1017 




1018 


81 


1019 




1020 




1021 




1022 




1023 


601 


1024 




1025 




1026 




1027 


802 


1028 




1029 




1030 


803 


1031 




1032 


805 


1033 




1034 




1035 


9009 


1036 


9000 


1037 


9001 


1038 




1039 


C 


1040 


85 


1041 


103 


1042 




1043 




1044 


C 


1045 


C 


1046 


86 


1047 


104 


1048 




1049 




1050 


C 


1051 


C 


1052 


87 


1053 


105 


1054 




1055 


C 


1056 


C 


1057 


1005 


1058 


1006 



WGH=( I (GHA+GHB)*.5)*DIFF)/TTIM 

GO TO 9001 

GHB*ELEV2 

WGH=( ( <GHAfGHB)*.5)*DlFF)/TTIM 

SGHsSGH+WGH 

GHA=ELEV2 

TIMZ=TIMEB 

RAD(4.END=107)NWELL.NCODENYRNMONDYNTIMiGHB 

IFINWELL-IWF.LLU 1)86.51.87 

CALL JULDY(NMONDY,NYRNTIM.TIMX) 

IF( T I MZ-T I MX 1201*20 1.200 

CALL JULDY(1231tNYR-1.2359.TM) 

DIFF=< I TM-TIMZJ+TIMX) 

GO TO 203 

DIFF=TIMX-TIMZ 

STIM=STIM+DIFF 

IF(STIM-TT IM)5580070 

WGH=I ( I (GHB+GHA)*.5)*DIFF)/TTIM) 

SGH=SGH+WGH 

GHA=GHB 

TIMZ=TIMX 

GO TO 50 

WGH=( < ( <GHB+GHA)*5)*DIFF)/TTIM) 

GHA=GHB 

TIMZ'TIMX 

GO TO 9001 

BACKSPACE 4 

BACKSPACE 4 

READI4,END=107)NWELL,NCODENYRiNMQNOY.NTIMELevl 

IFlNWELL-IWELLd M 86. 81 .8 7 

GHA=ELEV1 

CALL JULDY(NMONDYNYRNTIMtTIMEA) 

READ(4.ENDnl07)NWELL.NCODE.NYRlNMONDYNTIMELEV2 

CALL JULDYINMO.NDY.NYR1 NTIM*TIMEB) 

IFINYR1-NYR 1805 .805.801 

CALL JULDY(12.31NYR2359TM11 

DIFF=(TM1-T1MEA)+TIMEB 

SLOPF=(ELEV2-ELEV11/DIFF 

IF(JYR-NYR)803.803,802 

GHB=ELEV2-1SLOPE*(TIMEB-ETIME) ) 

WGH=( ( ( (GHB+GHA)*5HK(TM1-TIMEA) +ETIME) I/TTIM) 

GO TO 9001 

GHBELEV1+1SLOPE*(ETIME-TIMEAI ) 

GO TO 9000 

SLOPE=(ELEV2-ELEVn/(TIMEB-TIMEAl 

GHBa ELE V2- ( SLOPE* (TIMEB-E TIME) J 

IF (ETIME-TIMEA 19000. 9009* 9000 

ETIME=TIMEB 

WGH=( ( ( (GH3+GHA)*5t*lETIME-TIMF.A) J/TTIM) 

WGHP=SGH+WGH 

GO TO 1006 

WHEN NO RECORD IS FOUND FOR SELECTED TIME PERIOD THIS IS PRINTED IN OUTPUT 

WRITEI6tl03) I WELL! I ) t I MO i IDY t I YR . I TIMi JMOi JDY JYR JTIM 

FORMAT < ' ' ' WELL NOi ' .14.' NO RECORD FOR PERIOD OF ' 3I3i 15* ' TO 
1 ' 3I3. 15) 

GO TO 1005 

WHEN SELECTED TIME PERIOD STARTS BEFORE RECORDS ON TAPE THIS IS PRINTED IN 

OUTPUT* 

WRITEI6.1041 IWELH I ) 

FORMAT <' ( .' WELL NO. '>I4. 1 PERIOD OF RECORD STARTS BEFORE ACTUAL 
1RECORO' ) 

GO TO 1005 

WHEN SELECTED TIME PERIOD GOES BEYOND RECORDS ON TAPE THIS IS PRINTED IN 

OUTPUT 

WRlTE(6il05HWELL(I) 

FORMAT! ' '.' WELL NO. '.Ut* PERIOD OF RECORD BEYOND ACTUAL RECORD 
1' 1 

NUMBER 99999 STORED IN WATEL <KKK> INDICATES A MISSING RECORD FOR THESE 

SELECTED TIME PERIODS 

WGHP=99999. 

WATEL(KKK)WGHP 



1059 


1000 


JWELL1-IWELUI) 


1060 


1110 


CONTINUE 


1061 


107 


RETURN 


1062 




END 


1063 


c 




106* 


c 




1065 


c 


END OF RETRIEVAL PROGRAMS* STORAGE PROGRAMS FOLLOW. 


1066 


c 




1067 


c 


**** STORAGE PROGRAM FOR TEST HOLE DATA **** 


106* 




DEFINE FILE 1(21 ilOOOtU.IM) 


1069 




DIMFNSICN rtELLK 1000) (WELL 1 4 J NN 1 * J 


1070 




DIMENSION ISTART(M) 


1071 




DO 105 1=1.21 


1072 


105 


ISTART* I 1=0 


1073 


C 


READ FROM CARD THE NUMBER OF TEST HOLE VARIABLES TO BE STORED ON DISC 


107* 




fiEAO(5100! NOVAR 


1075 


100 


FORMAT* I*) 


1076 


c 


READ A VARIABLE NUMBER* FOUR VARIABLE VALUESi AMD FOUR TEST HOLE 


1077 


c 


MUMPERS FROM CARDS 


1078 


99 


READI5>200)NVAR i INN I I ) I 1 i* ) < WELL 1 I ) i 1 = 1 4 ) 


1079 


200 


FORMAT! l2*I5i*Fll0) 


1080 




IFlNVARU*il 


1031 


1 


IM-NVAR 


1082 




IF( ISTART(NVAR).EQ.O) GO TO 6 


1083 




REAOU'IMJ WELLl 


108* 




JM-JM-1 


1085 




GO TO 8 


1086 


6 


DO 7 I "it 1000 


1087 


7 


WELLlt I )0i 


1088 


8 


00 3 r*l* 


1089 




IFINNU ))33i2 


1090 


2 


IW=NN( I J 


1091 




WELLlt IWI-WELLI I 1 


1092 


3 


CONTINUE 


1093 


C 


WRITE TEST HOLE VARIABLE ON DISC 


109* 




WRITE 11 MM) WELL1 


1095 




rSTART(NVARJ=l 


1096 




GO TO 99 


1097 


C 


WRITE TEST HOLE VARIABLE DATA FROM DISC TO TAPE 


1098 


* 


IM-1 


1099 




00 5 I=liNOVAR 


1100 




READ11 ' IM) WELL1 


1101 




WRITE! 9) WELLl 


1102 


C 


WRITE VARIABLE NUMBER i TEST HOLE NUMBERi AND VARIABLE VALUE'ON PRINTER 


1103 




WRlTE1610m(JiWELLUJ)JlilOOO) 


no* 


101 


FORMAT t '0* 'VARIABLE NUMBER' tI2//>7U*2XFll*0)) 


1105 


5 


CONTJ.NUE 


1106 




WRITF(6102) 


1107 


102 


FORMAT Ml't'TEST HOLE VARIABLE DATA HAS BEEN STORED ON TAPE 1 ) 


1108 




STOP 


1109 




END 


1110 


C 




1111 


C , 




1112 


c 




1113 


c 


*** STORAGE PROGRAM FOR LAYER DATA **** 


11U 




DEFINE FILE 1 t 1000 , 161 U > IMJ 


1115 




DIMENSION STRATI 20, 91 tNWELLS ( 1000 ) 


1116 


c 


NOWEL IS THE NUMBER OF TEST HOLES FOR WHICH DATA IS STORED 


1117 




READ! 5 i 100) NOW EL 


1118 


100 


FORMAT* I*) 


1119 




ICOUN=0 


1120 




DO 3 IC*liNOWEL 


1121 


c 


NWELL IS THE TEST HOLE NUMBER 


1122 


c 


STRAT (1JI IS THE LAYER VARIABLES STORED IN EACH LAYER 


1123 


c 


MOLAY IS THE NUMBER OF THE LAYER FOR WHICH VARIABLES ARE STOREDi 


112* 




READ(5101)NWELLi(STRAT(l*Jt t Jl i 8 ) .NOLAY 


1125 


101 


FORMAT! Ut2F*.0lXiF5.0t2F3.0*2F7.0tF10i6I*J 


1126 




IF(WELL200i200i300 


1127 


300 


ICOUfi=ICOUfJ + l 


1128 




IFirtOLAY-U3.2l 


1129 


1 


RE ADI S 106) 1 (STRATI I *J) tJLi8) tI2tNOLAY) 



44 



mo 


106 


1131 


2 


1132 




1133 


3 


1134 


C 


1135 


200 


1136 




1137 




1138 




1139 


4 


1140 




1141 




1142 


5 


1143 


C 


1144 




1145 




1146 




1147 




1148 




1149 




1151 


103 


1152 


7 


1153 




1154 




1155 




1156 


104 


1157 




1158 




1159 


C 


1160 


C 


1161 


c 


1162 


c 


1163 


c 


1164 


c 


1165 


c 


1166 


c 


1167 


c 


1168 


c 


1169 


c 


1170 


c 


1171 


c 


1172 


c 


1173 


c 


1174 


c 


1175 




1176 




1177 




1178 




1179 




1180 




1181 




1182 




1183 




1184 




1183 


1 


1186 


c 


1187 


c 


1188 


c 


1189 


c 


1190 


5 


1191 


1000 


1192 




1193 




1194 




1195 




1196 




1197 


10 


1198 




1199 


c 


1200 


15 


1201 


c 



FORMAT{4X2F4.0ilX,F5.0t2F3.0t2F7.QiF10.6) 

WRITEIl'NWELDNOLAY.STRAT 

NWELLS< IO-NWEU 

CONTINUE 

SEQUENCE THE WELLS IN ASCENDING ORDER 

NUM=ICOUN-1 

DO 5 J=lt ICOUN 

DO 5 I=1NUM 

IFINWELLSf I >-NWELL5(I-t-l) >5i5t4 

HOLD=NWELLS< I ) 

NWELLSf NNWELLS1 1 + 1) 

NWELLSt I + 1UHOLD 

CONTINUE 

WRITF LAYER DATA ON TAPE AND PRINTER 

DO 7 1O1 ICOUN 

IM*NWELLSI 1C) 

READU' IM)NOLAYSTRAT 

IU-IM-1 

WRITE(S) lUtNOLAY.STRAT 

DO 6 I=liNOLAY 

FORMAT! ' ' t2I5i8F10l) 

CONTINUE 

IM*0 

WRITEUI IMiNOLAYtSTRAT 

WRIT(6104) 

FORMAT! 'I 1 t 'LAYER DATA STORED ON TAPE 1 ) 

STO.P 

END 



**## WATER LEVEL DATA STORAGE PROGRAM **** 
CALLS SUBROUTINE (TESTED 

PROGRAM STORES PUNCH CARD TAPE DOWN DATA AS WATER LEVEL ELEVATIONS ON 
MAGNETIC TAPE. 

MASTER TAPE OF WATER LEVEL ELEVATIONS IS CONSTRUCTED USING PIPE ELEVATIONS 
AND TAPE DOWN DATA STORED ON CARDS OR TAPES* 

WATER LEVEL DATA FOR SELECTED TEST HOLES CAN BE UPDATED USING THIS 
PROGRAM. 

TAPE PARAMETERS READ IN BY DATA STATEMENT 
WORDS-NUMBER OF WORDS WRITTEN PER RECORD 
LRECL-LENGTH OF RECORD SIZE 
BLKSIZ-BLOCK SUE PARAMETER 
BPI-BITS PER INCH FOR OUTPUT TAPE 
RGAP-RECORD GAPE BETWEEN RECORDS (INCHES) 
REAL*4 LRECLiINPBK 

DIMENSION KWL(M >ELVX(4I .ELVNI1000) i I TEST (1000) iJWELLO) 
DIMENSION IWN(100),IMO(100)iIDY!100)lYR(100) tELV(lOO) 
DIMENSION ITAPEOI iNMO ( 3 I .NDY (3 ) NYR ( 3 ) iNT IM( 3 } tEL i 3 ) 
COMMON NYRtNMOiNDYiNTIM.ELIYRtIMO.IDYfELVELVNJTESTtKTAPEt 
#NCODE I COUNT I PAGE I WN 

DATA WORDS>LRECL.BLKSIZtBPIRGAP/8..2<K2400. 1600.i.75/ 
DATA JWELL/3*0/ 
DO 1 I=l1000 
ITEST(I)1 
ELVNt I )"0.0 

READ IN TOP OF PIPE ELEVATIONS FOR ALL TEST HOLES. 
IV-LAST CARD TEST 
KWL m-TEST HOLE NUMBERS. 

ELVX (I)-TOP OF PIPE ELEVATIONS AT EACH TEST HOLE. 
READ(5tlOOOHV>(KWLU) illi4 t(ELVX( I) .I = l4) 
FORMAT! I2i4I5^Fll. 2) 
IF(IV.LE.O) GO TO 15 
DO 10 I=li4 

IF(KWLin.LE.O) GO TO 10 
KKK=KWL< I ) 
ELVN(KKK)ELVXU) 
CONTINUE 
GO TO 5 

READ IN PIPE ELEVATIONS CHANGES 
1*1 
IWN (I (-TEST HOLE NUMBER 



1202 


C 


IMO U)i IDY (I). IYR 11) IS THE MONTH. DAY, AND YEAR OF PIPE ELEVATION 


1203 


C 


IELV II)) CHANGES. 




1204 


16 


READ 15*1001) IWNI IJtlMOU I .IDYI I) >IYR(I ) ,FLVU) 




1205 


1001 


FORMAT<M5F10.2 ) 




1206 




IFI IWN< IULE.O) GO TO 17 




1207 




ITESTI IWNt I ) )0 




1208 




1 = 1 + 1 




1209 




GO TO 16 




1210 


17 


NOBS 1=1 




1211 


C 


NTAPE-NUMBER OF TAPE DRIVE UNITS USED. 




1212 


C 


ITAPE UI-TAPE DRIVE UNIT NUMBERS. 




1213 


C 


<TAPE-CREATED OR UPDATED TAPE UNIT NUMBER. 




121* 




READ I 5t 1002 ) M TAPE t II TAPE ( I ) . I=l,NTApE) .KTAPE 




1215 


1002 


FORMAT (51 5 




1216 


C 


READ WELL DUMBER FROM TAPES FOR I N I T I L I ZAT ION . 




1217 


1005 


00 1006 IT=1 iMTAPE 




1218 




ITPITAPEI IT) 




1219 




REAOI ITP. 103) JWELLI IT) 




1220 


1006 


CONTINUE 




1221 


1007 


WRITE 16*9999) 1 JWELLU ) , Iol,NTAPE) 




1222 


9999 


FORMAT! '0***** '3IlQt' **#**/) 




1223 


C 


FOLLOWING TEST STATEMENTS SELECT DATA FROM TWO OR THREE 


DATA SETS AND 


1224 


C 


ARRANGES THEM IN ORDER OF ASCENDING TEST HOLE NUMBERS FOR CREATING OR 


1225 


C 


UPDATING DATA TAPE. 




1226 




IFUWFLLU J .EQ0. AND JWELL(2) . EQ. 0. AND. JWFLL ( 3) .EO*0)GO 


TO 2000 


1227 




IF(JWFLLU) EQ.O. AND t JWFLL (2) tEQ.O.AND. JmFLL( 3) .GT.OIGO 


TO 1060 


1228 




IF<JWFLL(1).EQ.O.AND.JWFLL(2).GT.O.AND.JWFLL<3) .EQ.OlGO 


TO lObO 


1229 




!F( JWELL( D.GT.O -AND. JWELLI 2 J.COiO. AND. JWELLI 3) .EQ.OlGO 


TO 1030 


1230 




IF< JWELLU J .EQ.O. AND. JWELL(2) .EQ.JWELLI3) ) GO TO 1040 




1231 




IF(JWFLLUJ.EO.O.AND.JWELL(2).LT.JlvELL(3) GO TO 1050 




1232 




IF(JWELLU), EQ.O. AND. JWELLI2). GT.JWELLO) GO TO 1060 




1233 




!F(JV,'ELLI2), EQ.O. AND. JWCLLU).EQ.JWELL<3) GO TO 1020 




1234 




IFIJVELLI2). EG. 0. AND. JWELLU). LT.JWFLL(S) GO TO 1030 




1235 




IF<JWELLI2) .EQ.O. AND, JWELLU) ,GT. JWELLI 3) GO TO 1060 




1236 




ir<JWELL(3 .0.0. AND. JWELLll). EQ.JWELH2) GO TO 1010 




1237 




IFIJWF.LLI3 .0,0. AND. JWFLLID. LT.JWELL(2) GO TO 1030 




1238 




tFtJWELLO .EQ.O. AND, JWEILID.GT.JWELLI2) GO TO 1050 




1239 

1240 

1242 




IFIJWrLLll .EQ.JINELL<2I.AND.JWELL(2) EQ.JWELLOM GO TO 
IFUWELLI1 .EQ.JhELL<2).AMD.JWELL(2).LT.JWFLL<3)) GO TO 
FIJWzLLU -LT.JWELH2). AND. JWELLU). EQ.JWELLOI) GO TO 
IFIJWELLU -LT.JWELL(2).AND.JWELL(2).EQ.JWELL(3) ) GO TO 


1009 
1010 
1020 
1030 


1243 




IFIJWELLI1 -GT. JWELLU). AND. JWELL(2).EQ.JWELL<3M GO TO 


1040 


1244 
1245 
1246 
1247 
1248 
1249 




IFIJWELLU GT.JWELL(2).AND.JWELLI2).LT.JWFLL(3) ) GO TO 
IFIJWELLI1 .LT.JWELL<2).AND.JWLLI2).LT.JWELL(3) ) GO TO 
IFIJWFLLU .LT.JWELL<21.AND.JWEUL(2).GT.JWELL<3M GO TO 
IFIJWELLU GT.JWELL(2).ANf).JWELL(2).LT.JWELL(3M GO TO 

l I J R!rp^ L iiA^ TtJWELU(2t ' AND ' JVjFLL(21 * GT ' JhE L.L(3)) GO TO 
WK i I t ( o i 1 OOQ ) 


1050 
1030 
1030 
1050 
1060 


1250 
1251 


1008 


GS R TO T 2oSo >lTHERE ARE N TEST COMBINAT I^ FOR THESE TEST HOLES.'/) 


1252 


1009 


00 34 IT=ltNTAPE 




1253 




ITP=ITAPE(U) 




1254 




IFIITAPE(IT).EQ.O) GO TO 34 




1255 




BACKSPACE ITP 




1256 




NWFLL1=JWELLIITP) 




1257 


20 


[FUTP.EQ.U GO TO 21 




1258 
1259 


1 


^ A ^JJP' 103 'END=.30)NWELL ( NCODE. (NMOU ) .NDYII ) ,NYRU ) ,NT 


IM( I ) EL( I 


1260 




IREAD=3 




1261 


103 


FORMATtH,l2,3UX,3I2.I5tF7.2U 




1262 




GO TO 22 




1263 


21 


READ! I TP, 104 t END=30)NWFLI 




1264 


104 


FORMAT I l4t4I2l4*Ffl 71 'NCQDEtNMOU) NDY( 1) NYRI1) .NTI 


Mil) ,ELl 1) 


1265 




IREAD=1 * 




1266 


22 


IF(NWFLLUNE.NWELL) GO TO 25 




1267 
1268 




CALL TESTEL(NWELLiNOBSliIREAD) 




1269 




GO TO 20 




1270 


25 


JWELLI I TP)=NWELL 




1271 




GO TO 34 




.1272 


30 


JWRLLUTP)=0 
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1273 
1274 
1275 
1276 
1277 
1278 
1279 
1280 
12B1 
1282 
1233 
1284 
1285 
1286 
1287 
1288 
1289 
1290 
1291 
1292 
1293 
1294 
1295 
1296 
1297 
1298 
1299 
1300 
1301 
1302 
1303 
1304 
1305 
1306 
1307 
1308 
1309 
1310 
1311 
1312 
1313 
1314 
1315 
1316 
1317 
1318 
1319 
1320 
1321 
1322 
1323 
1324 
1325 
1326 
1327 
1328 
1329 
1330 
1331 
1332 
1333 
1334 
1335 
1336 
1337 
1338 
1339 
1340 
1341 
1342 
1343 



ITAPE1 IT) =0 

34 CONTINUE 
GO TO 1007 

1010 DO 45 IT=lt2 
ITP=ITAPE(IT) 

IF< ITAPE( IT).EQ0) GO TO 45 
BACKSPACE ITP 
NWELLl=JWcLL< I TP ) 

35 IF( ITP.EO.l) GO TO 36 

READl ITPtl03END=41)NWELLNCODEt ( NMO I I ) .NDYlI lNYR(I)iNTIM(l)ELII 
1) .1 = 1.3) 
IREAD=3 
GO TO 37 

36 READl ITP. 104 , END = 4l)NWELL , NCODEi NMOI 1) NDYU) NYR(1)NTIMU).EU1I 
IREAD=1 

37 IF(NWELLI.NE.NWELL) GO TO 40 
CALL TESTELINWELLiNOBSl.IREADI 



GO TO 35 
40 JWELL1 ITP)=NWtLL 

GO TO 45 
<tl JWFLLI ITP)=0 

ITAPE1 IT0 
45 CONTINUE 

GO TO 1007 
1020 DO 65 ITM.2 

GO TO (50.51) IT 

50 ITP=ITAPE<1) 

IFtlTAPEI 1) EQ.OIGO TO 65 
GO TO 52 

51 ITP=ITAPE<3) 
IF(ITAPEO). EQ.OIGO TO 65 

52 BACKSPACE ITP 
NWELLlaJWELLUTP) 

55 IFUTPiEQ.ll GO TO 56 

READi ITP103tEND*61)NWLL.NCODE (NMOt I 1 fNDYtl ) NYR(I ) .NT IM I I I iEL ( I 
U .1 = 1.3) 
IREAD=3 
GO TO 57 

56 READt ITP.104.END"6nNWELLNCODEiNM01 1 ) tNDY I 1 ) t.NYR ( 1 I .NT I M I 1 ) iEL ( 1 ) 
IREAD'l 

57 IF(NWFLLUNE.MWELL) GO TO 60 
CALL TESTEL1NWEUL.NOBS1.IREAD) 
NWELL1NWELL 

GO TO 55 

60 JWELLt ITP)=NWELL 
GO TO 65 

61 JWELL1 ITP)=0 
ITAPE1 ITP)=0 

65 CONTINUE 

GO TO 1007 
1030 ITPITAPEfll 

IF(ITAPEIl).EQ.O) GO TO 1007 
BACKSPACE ITP 
,NWELL1 = JWELL(ITP) 

70 IF(ITP.EQ.l) GO TO 71 

READl ITP.103END80)NWELL.NCOOEt t NMO ( I ) tNDY I I ) .NYRCI ) iNT IM( U t EL t I 
1) .1=1.3) 
IREAD=3 
GO TO 72 

71 READl ITP. 104, END = SO)NWELLNCODE. NMO U) .NDYUl NYR(U NT I M I 1 ) ,EL I 1 ) 
IREAD=1 

72 IFt.NWFLLltNE.NWELLI GO TO 75 
CALL TESTEHNWELLiNOBSliIREAD) 
NWELL1-NWELL 

GO TO 70 
75 JWELLdTPI-NWELL 

GO TO 1007 
80 JWELL1 ITP)*0 

ITAPE(1)=0 

GO TO 1007 
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1344 
1345 
1346 



1343 

1349 

1350 

1351 

1352 

1353 

1354 

1355 

1356 

1357 

1358 

1359 

1360 

1361 

1362 

1363 

1364 

1365 

1366 

1367 

1368 

1369 

1370 

1371 

1372 

1373 

1374 

1375 

1376 

1377 

1378 

1379 

1380 

1381 

1382 

1383 

1384 

1385 

1386 

1387 

1388 

1389 

1390 

1391 

1392 

1393 

1394 

1395 

1396 

1397 

1398 

1399 

1400 

1401 

1402 

1403 

1404 

1405 

1406 

1407 

1408 

1409 

1410 

1411 

1412 

1413 

1414 



1040 DO 100 ITlt2 

GO TOI85.86) ( IT 

85 ITPMTAPE<2) 
1FUTAPE(2).EQ0)GO TO 100 
GO TO 87 

86 ITPMTAPEO) 
IFUTAPEt31.EQ.OIGO TO 100 

87 BACKSPACE ITP 
NWELLl'JWELLUTP) 

90 IF(ITP.EQ.l) GO TO 91 
READ(ITP.103END=96)NWELLNCODE<NMOI 1 ) .NDYI I ) tNYR( I ) NT 1MU ) lEU I 

1).I=1,3) 
IREAD*3 
GO TO 92 

91 READ( ITP104iENO=96)NWELL.NCODENMO! 1 ) iNDY ( U iNYR( 1 ) iNTIMl 1 ) EU 1 ) 
IREAO=1 

92 IFtNKELLl.NE.NWELLl GO TO 95 
CAUL TESTEL(NWELLiNOSSl.lREAD) 
NhELLl-NWELL 

GO TO 90 

95 JKELLUTPI-NWELL 
GO TO 100 

96 JWELL(ITP)-0 
!TAPEIITP)0 

100 CONTINUE 

GO TO 1007 
1050 ITP-1TAPEI2J 

IF(ITAPE(2).EO.O)GO TO 115 
BACKSPACE ITP 
NWELLI-JWELLIITP) 

105 IFMTPiEQilt GO TO 106 

READUTPil03iENDlll)NWELLiNCODEi(NMO( I ) tNDY! I) NYR( I ) .NT IM ( It tEL ( 
Il)l*li3) 
IREAO=3 
GO TO 107 

106 READ(ITPil04iENDllllNWELL.NCODE.NMOIl) iNDY ( 1 ) NYR ( 1 ) .NT !M( 1) t EL 1 1 
*) 

IREAD=1 

107 IF(NWELLliNE.NWELL) GO TO 110 
CALL TESTELINWELL.NOBSltlREAD) 
NWELL1=NWELL 

GO TO 105 

110 JWfLL(ITPNWELL 
GO TO 115 

111 JWELL1ITP10 
ITAPE(ITP)=0 

115 CONTINUE 

GO TO 1007 
1060 ITP=ITAP(3) 

IFIITAPEOJ.EQ.OIGO TO 130- 
BACKSPACE ITP 
NWELL1=JWELL1ITP) 

120 IF(ITP.EO.l) GO TO 121 

READ(ITP.103iEND-126)NWFLLNCODE(NMO( I) .NOY ( I) tiNYR( I ) NT IM ( I) .EL ( 
II ) tl"l3) 
IREAO=3 
GO TO 122 

121 READUTPil04.ENO126)NWELLtNCODENMO(l) tNOYH) .NYR ( 1 I NT IM( 1 ) t EL I 1 

IREAD=1 

122 IFINWELL1.NE.NWELL) GO TO 125 
CALL TESTEUNWELLNOBSlilREAD) 
NWELL1=NWELL 

GO TO 120 

125 JWELLIITPI-NWELL 
GO TO 130 

126 JWELLIITP)0 
ITAPE(ITP)=0 

130 CONTINUE 

GO TO 1007 
: COMPUTES THE NUMBER OF FEET OF TAPE USED TC STORE WATER LEVEL ELEVATIONS 
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1415 
1416 
1417 
1418 
1419 
1420 
1421 
1422 
1423 
1424 
1425 
1426 
1427 
1428 
1429 
1430 
1431 
1432 
1433 
1434 
1435 
1436 
1437 
1438 
1439 
1440 
1441 
1442 
144-3 
1444 
1445 
1446 
1447 
1448 
1449 
1450 
1451 
1452 
1453 
1454 
1455 



C AND THE NUMBER OF RECORDS STORED. 
2000 RECD=ICOUNT 



INPBK=( (WORDS#BLKSU)/BPI I+RGAP 
FT( <RCD/RECDS)*INPBK)/12. 
WRITE(6i2001) I COUNT. KT APE FT 

2001 FORMAT! '0' * ' S$$S$ 'iI7' NUMBER OF OBSERVATIONS WRITTEN ON TAPE 
*' 5$iS$' /6Xi 'FEET USED TO STORE DATA ON THIS TAPE IS 'iF6.ll 
REWIND KTAPE 
STOP 
END 



SUBROUTINE TESTEL CALCULATES WATER LEVEL ELEVATIONS AND WRITES ON UPDATED 
TAPE. 

SUBROUTINE TESTELtNWELLiNOBSliIREADI 

COMMON NYR(3) NMO(3) iNDYO) tNTIM ( 3 ) EU I 3 ) * I YR I 100 ) 1 MO I 1 00 > I DYUO 
10) tELV(lOO) .ELVN(IOOO) i I TEST I 10001 .KTAPE .NCODE . I COUNT t I PAGE * IWNI 10 
10) 

DO 28 1=1 tIREAD 

IFdREAO.EQ.l ) GO TO 1000 

IFtNYRI I) .EQ.O) GO TO 28 

IF(ITEST(NWELL)*EQ.l) GO TO 27 

DO 25 Jl tNOBSl 

TEST FOR PIPE ELEVATION CHANGE 

GE. lYRUJ.AND.NMOl I I.GEtlMOf J) .AND.NDYI U .GE I DY ( J ) . AND . 



IFINYRI I ) 

1NWELL.EQ 

25 CONTINUE 

GO TO 27 

NEW PIPE 



IWNIJ) I GO TO 26 



ELEVATION 



26 ELVN1NWFLLI =ELVU) 
1TESTCNWELL1-1 

: WATER LEVEL ELEVATION CALCULATED. 

27 FL( I )=ELVNlNWELL)-EL( I ) 
OUTPUT WRITTEN ON UPDATED TAPE 

1000 WRITEIKTAPeilOAINhELLiNCODEiNMOt I 1 tNDYl I ) iNYRU I NT I Ml I ) EL I I) 
104 FORMAT! I4t4I2i I^F8.2 ) 

28 CONTINUE 
RETURN 
END 

: END APPENDIX C 
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APPENDIX D.-SOURCE LISTING FOR HYDRAULIC-COEFFICIENT 

RANGE-ASSIGNMENT PROGRAM 



1 


C 




2 


C 




3 


C 




4 


C 




5 


C 




6 


C 




7 


C 




9 


C 




10 






11 




200 


12 


C 


tt* 


13 






14 


C 


** 


15 






16 


C 


*# 


17 






IS 


C 


** 


19 






20 


C 


*# 


21 






22 




101 


23 


C 


## 


13 


C 


10 


26 






27 






26 


C 


** 


29 


C 


*# 


30 






31 






32 




20 


33 






34 






35 




30 


36 






3 / 






38 




40 


39 






40 






41 




50 


42 






43 






44 




500 


45 






46 


C 




47 


C 




48 


C 




49 




60 


50 






51 




70 


52 






53 




80 


54 






55 




90 


56 






57 




501 


sa 




100 


59 




110 


60 




1000 


61 




1002 


62 






63 




1001 


64 







IN EACH RANGE** 



APPENDIX D 



HYDRAULIC COEFFICIENT RANGE ASSIGNMENT PROGRAM 

***THIS PROGRAM ASSIGNS HYDRAULIC COEFFICIENT RANGE TO EACH LAYER*** 

THIS PROGRAM REQUIRES A SCRATCH TAPE FOR ASSIGNING CODE VALUES TO LAYER 

DATA 

DIMENSION AU35) ,A2(15) ,A3< 10) ,A4I15) >STRAT(208) A5I 15) 

**THE FOLLOWING CARD READS IN THE NUMBER OF MATERIAL TYPES 

READ15200)INOJNO,KNOLNOMNO 

FORMAT(5I5I 
(THE FOLLOWING CARD PUTS RANGE 1 MATERIAL INTO AN ARRAY FOR LOADING*** 

READ(5101) (Aid ) I = ltINOJ 
(THE FOLLOWING CARD PUTS RANGE 2 MATERIAL INTO AN ARRAY FOR LOADING*** 

REACH5.101) IA2(I).I = 1*JNO) 
(THE FOLLOWING CARD PUTS RANGE 3 MATERIAL INTO AN ARRAY FOR LOADING*** 

READ<5101) 1A3< I) I=1,KNO) 
(THE FOLLOWING CARD PUTS RANGE 4 MATERIAL INTO AN ARRAY FOR LOADING*** 

READiStlOl) (A4U )iI"lLNO 

MATERIAL 



'THE FOLLOWING CARD PUTS RANGE 
RFAD(5101)<A5(I)I=liMNO) 



INTO AN ARRAY FOR LOADING*** 



101 FORMATUOF8.0) 
***THE FOLLOWING CARD READS THE TEST HOLE NO., NO. OF LAYERS AND MATERIAL 



IN EACH RANGE*** 
STORED ON TAPE*** 



30,70,30 



40, 80*40 



TYPE 
REAOlBiEND=lQ01)NWELLtNOLAYiSTRAT 

DO 1000 I=l20 
IF (NOLAY.EG.O)GO TO 1002 

(THIS SERIES OF TESTS MATCHES MATERIAL TYPE TO THOSE 
WHEN MATERIAL TYPES MATCH THE RANGE IS ASSIGNED AND 
DO 20 J=liINO 

IF!STRAT(!,3)-A1<JM 20,60*20 
CONTINUE 
DO 30 J=liJNO 
IF(STRAT(Ii3)-A2UM 
CONTINUE 
DO 40 J=1,KNO 
IF(STRAT( Ii3)-A3(J) I 
CONTINUE 
DO 50 J=liLNO 

IF (STRATI I ,3)-A4(J) 150,90,50 
CONTINUE 
DO 500 J=liMNO 

IF I STRATI I ,3 )-A5<JM500 ,501,500 
CONTINUE 
GO TO 1000 

CODE 1.0, 2iO, 3tO, 4iO OR 5.0 IS ASSIGNED TO 
RANGE l.Oi 2.0, 3.0i 4.0, OR 5.0 RESPECTIVELY 
RANGE CRITERIA IS DETAILED EN TEXT 
CODE10 
GO TO 100 
CODE=2.0 
GO TO 100 
CODE=3*0 
GO TO 100 
CODE4.0 
GO TO 100 



DO 110 K=6>8 

STRATI I,K) = CODE 

CONTINUE 

WRITE(9)NWELL,NOLAY,STRAT 

GO TO 10 

REWIND 8 

ENDFILE 9 
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65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 



REWIND 9 

READ(9END=1003)NWELLNOLAYiSTRAT 
WRITE! B)NWELLNOLAY*STRAT 
GO TO 1004 
1003 REWIND 8 

WRITEI6. 10051 

STATEMENT INDICATES 
ON TAPE 1 ) 



C ***THE FOLLOWING 
1005 FORMAT ( ' ODATA 

STOP 

END 
C 



RANGES HAVE BEEN STORED ON TAPE*** 



END APPENDIX D 



APPENDIX E.-FLOW CHARTS AND LIST OF VARIABLES FOR 
PREPARING GRAPHICAL PRESENTATIONS 



Flow Charts 



MAP VIEW 



CROSS SECTION 



READ IN TEST HOLES TO BE 
PLOTTED 



READ SCALE FOR THE PLOT 

XSCAL.YSCAU 

(SEE VARIABLE LIST) 



READ SPECIFIED BOUNDARY 
OF MAPPED AREA FROM 
DISK FILE 



SCALE FIELD DATA POINTS 
TO PLOTTER COORDINATES 
(SYSTEM DERNED BY USERS! 



PLOT MAP OUTLINE 



SUBROUTINE WELPT 
READ NW, ALON, ALAT 
FROM CARDS OR DISK RLE 
[SEE VARIABLE LIST) 



PLOT SELECTED TEST 
HOLES 




SCALE PARAMETERS 
XLENG, YLENG , GTEST, BTEST 
(SEE WRIA8LE LIST) 



RETRIEVE LAYER \ttHIABLES 

AND STORE ON DISK , NWELL, 

THKK.KVWJJ 

(SEE TRIABLE LIST) 



RETRIEVE TEST-HOLE 
WBABLES 8 STORE ON DISK 
NWELL, OSURF, TDEPH, MAXU 
(SEE VARIABLE LIST) 



CALL SCALE PLOT 
SUBROUTINE 




CALL GRID AXIS SUBROUTINE 
AND ANNOTATE THE ORID 



PLOT INDIVIDUAL TEST-HOLE 
LOO AND WATER SURFACE 



FIGURE E-l.-Flow charts for preparing map views and cross sections. 
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List of Variables 



Variable 
name 



XLENG 
YLENG 
GTEST 

BTEST 



NWELL 
IVAR 

THICK 
KVALU 



NWELL 
GSURF 



TDEPH 

MAXLA 



Format Deso iption of variable 
Boundary Card for Cross Section 

F10.0 Distance in inches between 
test-hole-log plots. 

F10.0 Length in inches of the cross- 
section plot. 

F10.0 Reference ground-surface 

datum from which logs are 
measured on the plot. 

F10.0 Lower reference bedrock 

datum beneath which test- 
hole logs do not extend. 



Lithology Card for Cross Section 

14 Test-hole numbers of logs 

used in cross section, 
14 Layer number for lithology 
of test holes in cross sec- 
tions. 

F9.1 Thickness of lithologic layers 
in feet. 

F9.1 Hydraulic-coefficient range 
(1-5) of the lithology in 
each layer. 

Test-Hole Card for Cross Section 

18 Test-hole number for each 
log in cross section. 

F12.3 Ground-surface elevation 

(feet above mean sea 
level) for the top of each 
log. 

F12.3 Total depth of each test hole. 
12 Number of lithology layers in 
each test hole in cross sec- 
tion. 



Variable 
name 



Format Description of vai iable 
Location Card for Map Views 



Ground-Water-Elevation Card for Cross Section 

HWELL 15 Test-hole number. 

HCODE 12 Code indicating the method 

of measurement of ground- 
water level (1 continuous 
recorder, 2= manual tape 
down) , 

NYR 13 Year when water surface was 

recorded. 

NMO 13 Month when water surface 

was recorded. 

NDY 13 Day when water surface was 

recorded. 

NTIM 15 Time (military) when water 

surface was recorded. 

ELEV F8.2 Ground-water-surface eleva- 

tion (feet above mean sea 
level) in the test hole. One 
card per test hole. 



Note: four maps are plotted: distributive map, sub- 
datum map, hydraulic-coefficient map, isopachous 
map. 

NW F12.3 Test-hole number for each 

hole plotted in map view. 

ALON F12.3 Longitude in degrees, min- 

utes, seconds of test-hole 
locations. 

ALAT F12.3 Latitude in degrees, minutes, 

seconds of test-hole loca- 
tions. 

Lithology Card for Subdatum Map 

NW 13 Test-hole number of each 

test hole that penetrates 
the specific datum plane. 

ISUB 16 Elevation of the datum plane 

to be plotted (feet above 
mean sea level). 

IDEPH 14 Depth below ground surface 

to the top of the lithologic 
layer in which specified 
datum plane lies. 

ELTOP F8.2 Elevation of the top of the 

lithology layer that con- 
tains the specified datum 
plane (feet above mean sea 
level). 

RANGE F12.4 Range of hydraulic coeffi- 

cients (1-5) at the speci- 
fied datum. 

THICK F12.4 Thickness in feet of the lith- 

ology layer found at the 
specified datum. 

DAPLT F12.4 This variable name is as- 

signed to the range of K, 
S t or T to be plotted. 

Test-Hole Card for Map View 

Note: four maps are plotted: distributive map, sub- 
datum map, hydraulic-coefficient map, isopachous 
map. 

NW 16 Test-hole number. 

THVAR F5.0 Variable THVAB may be 

any of the 21 test-hole 
variables used for plotting. 

Net-Thickness Card for Isopachous Map 

WELL NO. A10 Descriptive literal field. 

NW 16 Test-hole number. 

NET THICK A10 Descriptive literal field. 

THICK F13.0 Net thickness in feet of 

lithologic layers having the 
selected hydraulic-coeffi- 
cient ranges. 
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ana e Format Description of vanable "''" 6 Format Description of variable 



)lflttlC 

Net-thicknes,s Card for Isopachous Map Continued Card for Hydraulic-Coefficient Map Continued 

CODES AC Descriptive literal field. AKVAL F5.0 Permeability averaged for all 

layers in the test hole, gpd/ 
CODE (I) 5F4.0 Number (1-5) of the ranges ftj flf aquiferi The slzo of 

of hydraulic coefficients, the fields containing liy- 

thal is, K, S, or T, to be draulic coefficients may be 

plotted. alteied to accommodate re- 

trieved data. 
ASVAL F5.0 Storage coefficient, dimen- 

Card for Hydraulic-Coefficient Map sionless value. 

ATVAL Fll.O Transmissibility, gpd /ft of 

NW 14 Test-hole number, aquifer. 
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APPENDIX F.-INPUT CARDS FOR GRAPHICAL PRESENTATIONS 



Stj-tifi 
V1XW 



9E-S2 
HdBQl 



ddnss 



8-t 

TI3HN 




tt-ll 



\Z-tt 
HUN 



9T-bI 
AQN 



ET-IT 
OWN 



OI-B 
BMI 



3QODN 



S-T 
TI3HN 



^jgs 

iSu,S 



1S31B 
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504 S 70 514 1346 1020.10 
AB C E F 



GROUND-WATER- 
ELEVATION CARD 



504.000 1037.300 43.000 

ABC 



S.OOO 



TEST-HOLE CARD 



BLANK CARD 



604 I 8.00 2.00 

AS C D 



LITHOLOGY CARD 



20. 
B 



1080. 
C 



940. 




BOUNDARY CARD 



FIGURE F-3 Input-card setup for plotting cross sections. 



504 43 
A B 



TEST-HOLE DATA CARD 



504000 974839.000 345751000 
ABC 



LOCATION CARD 



FIGURE F-4. Input-card setup for plotting test-hole data (see figure 10} . 
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LITHOLOGY DATA CARD 





519.000 981205000 350327.000 
ABC 



LOCATION CARD 



FIGURE F-6. Input-card setup for plotting data at a specified subdalum 

(see figure 12). 



504 293 003 
ABC 



HYDRAULIC- 
COEFFICIENT CARD 



504 000 974839.000 345751.000 
ABC 



LOCATION CARD 



FIGURE F-G. Input-card setup for plotting average peimeability coefficients 

(see figure 14). 



WELL NO. 
A 


S04NET 
B 
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C 
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D 
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H 1 J 
THICKNESS 


CARD 





504.000 974639,000 345761.000 
ABC 



LOCATION CARD 



FIGURE F-7. Input card setup for isopachous maps (see figure 16) . 
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TRADE NAMES ARE USED IN THIS PUBLICATION SOLELY FOR THE PURPOSE OF PROVIDING SPE- 
CEFIC INFORMATION. MENTION OF A TRADE NAME DOES NOT CONSTITUTE A GUARANTEE OR 

WARRANTY OF THE PRODUCT BY THE u.s DEPARTMENT OF AGRICULTURE OR AN ENDORSE- 
MENT BY THE DEPARTMENT OVER OTHER PRODUCTS NOT MENTIONED 
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